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This is an appeal from an Office Action dated March 9, 2005 in which claims 1, 2 and 5 
were finally rejected. 

RK AT. PARTY TN INTEREST 

Seagate Technology LLC, a limited liability corporation organized under the laws of the 
state of Delaware, and having offices at 920 Disc Drive, Scotts Valley, California 95066, has 
acquired the entire right, title and interest in and to the invention, the application, and any and all 
patents to be obtained therefor, as set forth in the Assignment filed with the patent application and 
recorded on Reel 01 1063, frame 0119. 
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STATTTS OF THF. CLAIMS 

L Total number of claims in the application. 



Claims in the application are: 8 
n. Status of all the claims. 

A. Claims canceled: 0 

B. Claims withdrawn but not canceled: 0 

C. Claims pending: 8 

D. Claims allowed: 0 

E. Claims rejected: 3 

F. Claims objected to: 5 
in. Claims on appeal 

The claims on appeal are: 1-8 



STATTfS OF AMKNDMF.NTS 

An Amendment was filed after the mailing of the final rejection. The Amendment was in 
response to the 35 USC § 1 12, first paragraph and 35 USC § 1 12, second paragraph rejections. The 
Examiner indicated, in the Advisory Action, that the Amendment was entered. Therefore, the 
Appellant believes the § 112 rejection was withdrawn. 

STIMMARY OF CLATMFft STIBtTFd MATTER 
L INTRODUCTION 

The present invention relates to a disc drive. More particularly, the present invention relates 
to a process of designing an optimal passive damping mount for a disc drive. 

II. BRIEF BACKGROUND 

There are two types of disturbances that affect operation of disc drives. These two types of 
disturbances are internal disturbances and extemal disturbances. Examples of intemal disturbances 
include wdndage created by the rotating disc and encountered by the head/arm assembly, noise in 
the position error signal (PES), structural resonance of the disc drive, motion of the actuator 
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assembly, imbalance of the discs, written-in run out in the servo patterns due to disturbances during 
writing of the servo patterns, unmodeled dynamics and non-linearities, such as pivot bearing 
friction. Examples of external disturbances include environmental (seismic) vibration and external 
reactions to intemal disturbances. 

When a disc drive is subjected to disturbances, head/disc tracking is adversely affected. 
Generally, the heads and the data tracks on the disc have some defined width that is large enough to 
assure that disturbances do not materially affect the performance of the disc drive. However, the 
defined width of the data tracks limit the maximum areal data density of the disc. Typically, passive 
damping mounts are mounted to the bottom of the disc drive to isolate disturbances. As the stif&iess 
of the damping mounts decreases, better attenuation of seismic vibration is obtained. However, soft 
isolation mounts cannot attenuate internally generated disturbances. Generally, internally generated 
disturbances can be attenuated by hard or stiff mounts. Therefore, external and intemal vibrations 
have conflicting effects when using passive damping mounts. Attenuation of one source of 
disturbance by passive mounts leads to amplification of the other source. 

IIL THE PRESENT INVENTION 

A. INDEPENDENT CLAIM 1 

Independent claim 1 of the present invention provides a process of designing an optimal 
vibration mount for a disc drive that eliminates the conflicting effects of a mount that is too stiff or 
a mount that is too soft. The process includes the steps of computing (202) an extemal disturbance 
model for the disc drive and computing (202) an intemal disturbance model for the disc drive (see 
the Specification on page 5, line 25 through page 7, line 4 and page 9, line 15 through line 19). The 
process includes defining (200) an inertia matrix for the disc drive (see Specification on page 9, 
lines 15 through 16). The process includes defming (204) a state estimator based on the inertia 
matrix and extemal and intemal disturbance models to minimize a defined norm of a state 
estimation error (see the Specification on page 7, line 5 through page 8, line 4 and page 9, line 20- 
25). The process includes calculating (206, 208) the gain of the state estimator as a solution to a 
filter algebraic Riccati equation (see the Specification on page 8, line 5 through page 9, line 12 and 
page 10, line 1 through line 6). The process also includes defining (210) optimal mount damping 



and stiffiiess parameters based on the calculated state estimator gain (see the Specification on page 

10, lines 6 through 9). 

CiROTINDS OF RKTF.CTTON TO BE RKVTEWF.D ON APPKAT . 

L Whether claims 1-2 are unpatentable over Kim et al. (See Appendix B, Exhibit A) 
in view of Brown et al. (See Appendix B, Exhibit B) or Lee (See Appendix B, Exhibit C). 

n. Whether claim 5 is unpatentable over Kim et al. in view of Brown et al. and further 
in view of Reif et al. (See Appendix B, Exhibit D) or Shah et al. (See Appendix B, Exhibit E). 

ARr.TTMFNT 

L INTRODUCTION 

The present invention provides a process of designing an optimal vibration mount or 
optimal passive damping mount for a disc drive. Typically vibration mounts or damping mounts are 
used to attenuate seismic or external disturbances as well as to attenuate internal disturbances of a 
disc drive. Soft mounts can attenuate seismic disturbances well, while stiff mounts can attenuate 
internal disturbances well. These conflicting effects of soft mounts versus stiff mounts causes 
attenuation of one source of disturbance and amplification of the other source of disturbance. 
Embodiments of the present invention include steps for designing an optimal mount that can 
attenuate both intemal and extemal disturbances without producing conflicting effects. Claims 1, 2 
and 5 are currentiy rejected by at least three different references. 

11. PROSECUTION OF THE PRESENT APPLICATION 

Appellant filed patent application serial number 09/651,046 on August 30, 2000, claiming 
priority to provisional application serial nximber 60/162,269 filed October 28, 1999. The original 
non-provisional application set forth claims 1-8, which are still pending. In an Office Action dated 
July 12, 2004, claims 1, 2 and 5 were rejected and claims 3-4 and 6-8 were objected to. Appellant 
responded on October 12, 2004 by way of an Amendment that amended claim 1. Thereafter, 
Appellant received a Final Office Action mailed March 9, 2005 rejecting claims 1-8 under 35 USC 
§ 1 12, first and second paragraphs, rejecting claims 1, 2 and 5 under 35 USC § 103(a) and objecting 
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to claims 3-4 and 6-8. Appellant responded to the Final Office Action on May 5, 2005 by way of an 
Amendment that amended claim 1. The Amendment After Final amended claim 1 back to its 
original form. The Examiner indicated in the Advisory Action mailed June 3, 2005 that the 
Amendment After Final was entered. Therefore, the 35 USC § 1 12 rejection should be considered 
withdrawn. However, claims 1, 2 and 5 still remain rejected under 35 USC § 103(a) and claims 4-6 
and 7-8 still remain objected to. 

III. PRIOR ART 

Kim et al. (Exhibit A) teaches a control system design method that uses Linear Matrix 
hiequality (LMI), that replaces or complements the Lyapunov-Riccati equations, to suppress shaft 
vibration of an industrial motor drive system. The resulting method proposes an H infinity speed 
controller that uses weights to satisfy tracking and disturbance rejection for an induction motor. The 
control system design of the Kim et al reference includes the use of three different concepts. These 
are: a two-mass resonant system, an H infinity generalized plant and a control system design by 
LMI. The Kim et al. reference validates the proposed controller by comparing it to a conventional 
state feedback controller using a pole placement method. The Examiner further attached to the Kim 
et al. reference an Abstract of a different paper that the same author had published. However, the 
Examiner neither cited this as a reference in the PTO-892 form nor referred to it in any of the 
mailed Office Actions as a reference. Therefore, although the Appellant acknowledges the 
attachment to the Kim et al. references, the Appellant is not addressing its existence in the 
following arguments. 

Brown et al. (Exhibit B) teaches the design of a research servowriter to record magnetic 
position encoding schemes in magnetic disk storage applications and to evaluate disk drive 
technologies. The servowriter includes signal processing hardware, a programmable "pattern 
generator" controlled by a supervisory computer and a laser interferometric system to measure the 
position of the read/write head position to thereby estimate the velocity of the read/write head for 
implementation of a state-variable-based proportional-integral-derivative controller on a 
microprocessor. The reference also discusses error sources. Sources of error include environmental 
factors, such as seismic vibration and humidity, spindle and head suspension jitter and other 
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external vibrations. The reference suggests selecting in-line heads to reduce tangential jitter, 
isolating the servowriter from external vibrations by mounting it on a large granite base that floats 
on air isolation mounts and placing the servowriter in a clean room to prevent environmental 
contamination. 

Lee (Exhibit C) proposes an H infinity model matching algorithm for motor tracking 
control. The control method of Lee provides improved drive performance and reduced sensitivity of 
disturbances. The algorithm is based on the discrete Riccati equation and the H infinity 
interpolation theory. The controller effectively controls the presence of external disturbances. 

Reif et al. (Exhibit D) teaches non-linear state observation using a Riccati differential 
equation to improve the extended Kalman filter. The Reif et al reference applies the non-linear state 
observation to an induction motor for the estimation of flux and angular velocity, hi addition, the 
Reif et al. reference discusses the proposed non-linear state observation in terms of the robust 
control theory and H infinity filtering. 

Shah et al. (Exhibit E) teaches a method for estimating state variables of non-linear control 
systems using improved extended Kalman filtering techniques. The improved Kalman filtering 
techniques include using a continuous-time extended Kalman filter (CTEKF) to estimate all state 
variables. The CTEKF uses discrete time measurement and integration. 

IV. THE REJECTION OF CLAIMS 1, 2 AND 5 SHOULD BE REVERSED 

Claims 1, 2 and 5 have been rejected under 35 USC § 103(a) with various combinations of 
the cited references. To establish a prima facie case of obviousness, three basic criteria must be 
met. First, there must be some suggestion or motivation either in the references themselves or in 
the knowledge generally available to one of ordinary skill in the art, to modify the reference or to 
combine reference teachings. Second, there must be a reasonable expectation of success. Finally, 
the prior art reference (or references when combined) rmi<;t teach nr <;iiggest all of the claim 
limita ti ons. Tn re V a ec k, 20 U.S.P.Q.2d 1438 (Fed. Cir. 1991); M.P.E.P. §2143. 

Li the Response to Argimients section of the Final Rejection, the Examiner stated that 
"Applicants have merely paraphrased selected portions of the applied art without addressing the 
sections of the applied art as indicated in the last Office Action." It is respectfully pointed out 
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that in accordance with MPEP § 2142, "The examiner bears the initial burden of factually 
supporting any prima facie case of obviousness. If the examiner does not produce a prima facie 
case, the applicant is under no obligation to submit evidence of non-obviousness." It is 
respectfully submitted that where no prima facie case of obviousness is substantiated, the 
Applicant has no burden to reply. However, the Appellant, as will be thoroughly discussed 
below, has diligently pointed out elements of the claims that are not taught or suggested by the 
cited references. Furthermore, 37 CFR §1.1 04(c)(2) states the following: 

In rejecting claims for want of novelty or for obviousness, the examiner must cite the best 
references at his or her command. When a reference is complex or shows or describes 
inventions other than that claimed by the applicant, the particular art relied on must be 
designated as nearly as practical. The pertinence of each reference, if not apparent, must 
be clearly explained and each rejected claim specified. 

Not only has the Examiner not satisfied the burden of showing all of the elements of the claims 
in the cited references, the Examiner has cited references which have features that are not readily 
apparent as being relevant to elements of the claims. The cited references are complex and are 
not pertinent to the scope of the present invention. 

A- riflim^; 1 and ? are Allowahle Over the Art of Record 

Sections eight and nine of the Final Rejection indicated that claims 1 and 2 were rejected 
under 35 USC § 103(a) as being xmpatentable over Kim et al. in view of Brown et al. or Lee. 
Appellant respectfully submits that the evidence of record in the Final Rejection is insufficient to 
establish a prima facie case of obviousness. In particular, it is respectfully submitted that even if the 
Kim et al. and Brown et al. or Lee references were combined, the combination of cited references 
simply do not teach or suggest all of the elements recited in claims 1 and 2. In particular, the 
combination of references fail to teach or suggest: 

■ computing an external disturbance model for the disc drive 

■ computing an internal disturbance model for the disc drive 

■ defining an inertia matrix for the disc drive 
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■ defining a state estimator based on the inertia matrix and external and internal 
disturbance models to minimize a defined norm of a state estimation error 

■ calculating the gain of the state estimator as a solution to a filter algebraic Riccati 
equation 

■ defining optimal mount damping and stifl&iess parameters based on the calculated state 
estimator gain 

Before thoroughly discussing each of the above bulleted elements, it should be noted that 
the rotation drives and drive systems that are referred to in the Kim et al. and Lee references are 
drives that rotate and include means for giving motion to a machine or part of a machine. The 
rotation drives and drive systems referred to are not a disc drive as claimed in claim 1 . As described 
in the Specification with reference to FIG. 1, a disc drive is a device for storing data. Although a 
spindle motor is located in a disc drive for rotating a disc, a disc drive includes many other 
components for reading and vmting data to a disc. These components and their operations are 
discussed in connection with FIG. 1 of the Specification. 

In section nine of the Final Rejection, the Examiner fmds that Kim et al. discloses 
"computing an external disturbance model for a disc drive" and "computing an internal disturbance 
model for a disc drive" as recited in independent claim 1 . It is respectfully submitted that the Kim et 
al. reference neither discloses modeling external disturbances for a disc drive nor modeling intemal 
disturbances for a disc drive. In fact, the Kim et al. reference does not discuss extemal disturbances 
at all and, instead, teaches a method of controlling vibration of a drive shaft. Instead, and as 
mentioned above, the Kim et al. reference teaches a method of designing a control system for a 
motor drive system that suppresses shaft vibration. The Specification of the present invention 
specifically notes that examples of extemal disturbances for a disc drive are seismic vibrations that 
are produced extemal to the disc drive. Motor shaft vibration is not an extemal disturbance to a 
drive system. In addition, the method taught in the Kim et al. reference is specific to vibration 
suppression in an industrial or induction type motor. Kim et al. does not teach or suggest computing 
a model for any type of disturbance (externally or intemally) for a disc drive. In fact, the Kim et al. 
does not even mention a disc drive. Although shaft vibration can potentially be a source of intemal 
disturbance in a disc drive because a disc drive includes a spindle motor, it is well knovra, and 
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discussed in the Background of the Specification, that disc drives include many more sources of 
internal disturbances that are not considered in the Kim et al. reference. 

Furthermore, neither the Brown et al. reference nor the Lee reference correct deficiencies of 
the Kim et al reference as pertaining to computing "an external disturbance model for a disc drive" 
and "an intemal disturbance model for a disc drive". Brovm et al. does not teach or suggest 
modeling extemal disturbances for a disc drive nor modeling intemal disturbances for a disc drive. 
Although Brown et al. identifies error sources for a servo writer. Brown et al. suggests ways of 
mitigating the identified error sources by selecting in-line heads to reduce tangential jitter, isolating 
the servowriter from extemal vibrations by mounting in on a base that floats on air isolation mounts 
and placing the servowriter in a clean room. Brown et al. does not attempt to model intemally or 
externally derived disturbances for a disc drive. In addition, the Lee reference does not teach or 
suggest modeling extemal or intemal disturbances for a disc drive. Although Lee identifies load 
torque on the motor shaft of a motor drive system as being a disturbance to consider, Lee does not 
disclose modeling disturbances for a disc drive. 

In section nine of the Final Rejection, the Examiner finds that Kim et al. defines "an inertia 
matrix for the disc drive" as recited in independent claim 1 . However, It is respectfully submitted 
that Kim et al. does not teach or suggest such a claim element. Kim et al. discloses a motor inertia 
value and a load inertia value when providing a two-mass resonant system. Although a disc drive 
includes a spindle motor, an inertia matrix for a disc drive, in its entirety, is not defined in Kim et 
al. As discussed in the Specification in reference to FIG. 1, a disc drive includes many more 
components than just a motor. Such components are not considered by the Kim et al. reference. 

Furthermore, neither the Brown et al. reference nor the Lee reference correct deficiencies of 
the Kim et al reference as pertaining to "defining an inertia matrix for the disc drive". Neither 
Brown et al. nor Lee teach or suggest defining an inertia matrix for a disc drive. 

In section nine of the Final Rejection, the Examiner finds that Kim et al. defines a "state 
estimator based on the inertia matrix and extemal and intemal disturbance models" and calculates 
"the gain of the state estimator as a solution filter algebraic Riccatti equation" as recited in claim 1 . 
It is respectfully submitted that the Kim et al. reference, as discussed above, does not teach or 
suggest such claim elements. It is believed that the Kim et al. reference does not teach or suggest 
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defining an inertia matrix or computing extemal and intemal disturbance models for a disc drive. 
Therefore, Kim et al. could not possibly define a state estimator based on these elements as well as 
calculating a gain of such a state estimator. The Examiner states that "Kim et al. discloses use of 
Riccati solutions to the reduction of vibration in rotational drives including inertial matrices for 
intemal disturbances" in section 2-4 of the reference. This statement and the information relayed in 
sections 2-4 of the reference do not read on these or any of the above listed claim elements. 

Furthermore, neither the Brown et al. reference nor the Lee reference correct deficiencies of 
the Kim et al reference as pertaining to "defining a state estimator based on the inertia matrix and 
extemal and intemal disturbance models to minimize a defmed norm of a state estimation error" 
and "calculating the gain of the state estimator as a solution to a filter algebraic Riccati equation". 
Neither Brown et al. nor Lee teach or suggest defining a state estimator based on the inertia matrix 
and extemal and intemal disturbance models to minimize a defined nomi of a state estimation error 
or calculating the gain of the state estimator as a solution to a filter algebraic Riccati equation. 

In section nine of the Final Rejection, the Examiner finds that Kim et al. defines "optimal 
mount damping stif&iess parameters based on the calculated state estimator gain" as recited in 
claim 1. It is respectfiiUy submitted that the Kim et al. reference neither discloses defining optimal 
mount damping stifi&iess parameters nor defining optimal mount damping stiffiiess parameters 
based on the calculated state estimator gain. In fact, the Kim et al. reference does not discuss 
mounts for a disc drive at all and, instead, teaches a control system for suppressing vibration of a 
motor shaft. In addition, it is believed that the Kim et al. reference does not teach or suggest 
defining an inertia matrix, computing extemal and intemal disturbance models for a disc drive, 
defining a state estimator, or calculating the gain of the state estimator. Therefore, Kim et al. could 
not possibly calculate any parameter, let alone an optimal mount damping and stifi&iess parameter, 
based on the state estimator. 

Furthermore, neither the Brown et al. reference nor the Lee reference correct deficiencies of 
the Kim et al reference as pertaining to "defining optimal mount damping and stiffiiess parameters 
based on the calculated state estimator gain". Although Brown et al. identifies the use of a base that 
floats on air isolation mounts. Brown et al. does not teach or suggest a method of defining optimal 
mount damping and stiffiiess parameters nor teaching or suggesting a method of defining optimal 



mount damping and stiffiiess parameters based on the calculated state estimator gain. Lee does not 
teach or suggest a disc drive, let alone optimal damping mounts. 

In section ten, the Final Rejection conceded that Kim et al. did not disclose taking into 
account external disturbances in addition to the intemal disturbances. The Final Rejection turns to 
the teachings of Brown et al. or Lee as disclosing extemal as well as intemal disturbances to motor 
drives. In particular, the Final Rejection foimd that Brown et al. identifies environmental factors 
such as seismic vibration, temperature change and steps to contain these error sources are described. 
The Appellant finds this concession irrelevant and moot in view of the fact that the Examiner has 
not provided a prima facie case of non-obviousness. Specifically, the Examiner has not shown that 
other elements of the claim (as bulleted above) are found within the combination of references. 

In section nine of the Final Rejection, the Examiner finds that Kim et al. calculates 
"covariance matrix based on the solution to the filter algebraic Riccati equation" and calculates "the 
state estimator gain based on the covariance matrix" as recited in claim 2. In one aspect, it is 
respectfiilly submitted that Kim et al. does not teach or suggest the elements of claim 2. Kim et al. 
does not disclose the use of a covariance matrix, let alone the use of a covariance matrix to 
calculate the state estimator gain. In another aspect, claim 2 is allowable over the cited reference as 
depending fi-om an allowable base claim. 

Furthermore, neither the Brown et al. reference nor the Lee reference correct deficiencies of 
the Kim et al reference with respect to the elements of claim 2. Neither Brown et al. nor Lee teach 
or suggest calculating a covariance matrix based on the solution to the filter algebraic Riccati 
equation nor calculating the state estimator gain based on the covariance matrix. 

It is respectfiilly believed that claims 1 and 2 are in condition for allowance. Reversal is 
respectfiilly requested. 

K Claim 5 is Allowable Over the Art of Record, 

In section fifl:een of the Office Action, claim 5 is rejected under 35 USC § 103(a) as being 
unpatentable over Kim et al. in view of Brown et al. and fiirther in view of Reif et al. or Shah et al. 
Appellant respectfiilly submits that the evidence of record in the Final Rejection is insufficient to 
establish a prima facie case of obviousness. In particular, it is respectfiilly submitted that even if the 
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Kim et al. and Brown et al. and further the Reif et al or Shah et al. references were combined, the 
combination of cited references simply do not teach or suggest all of the elements recited in claim 
5. In particular, the combination of references fail to teach or suggest that the state estimator of 
claim 1 is a Kalman Filter as claimed in claim 5. 

The Examiner concedes, in section sixteen of the Final Rejection, that Kim et al. is silent on 
whether the state estimator is a Kahnan Filter. In section seventeen, the Examiner looks to Reif et 
al. or Shah et al. as teaching a Kalman Filter. The Appellant finds this concession irrelevant and 
moot in view of the fact that the Examiner has not provided a prima facie case of non-obviousness 
related to other elements included in claim 5 and as discussed above with respect to claim 1 . In 
particular, claim 5 is allowable over the cited references as depending on an allowable base claim. 

V. CONCLUSION 

The cited references in the Final Rejection fail to teach or suggest the feature recited in 
claims 1, 2 and 5. Accordingly, Apellant respectfully requests that the Board reverse the Examiner 
and find that claims 1, 2 and 5 are in condition for allowance. 



Respectfully submitted. 



WESTMAN, CHAMPLIN & KELLY, P.A. 




Suite 1400 - International Centre 
900 Second Avenue South 
MinneapoUs, Minnesota 55402-3319 
Phone: (612)334-3222 
Fax: (612)334-3312 
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Appendix A: Claims 

A process of designing an optimal vibration moimt for a disc drive comprising steps of: 

a) computing an external disturbance model for the disc drive; 

b) computing an internal disturbance model for the disc drive; 

c) defining an inertia matrix for the disc drive; 

d) defining a state estimator based on the inertia matrix and external and intemal 

disturbance models to minimize a defined norm of a state estimation error; 

e) calculating the gain of the state estimator as a solution to a filter algebraic Riccati 

equation; and 

f) defining optimal mount damping and stiffness parameters based on the calculated 

state estimator gain. 

The process of claim 1, wherein step (e) includes, 

el) calculating a covariance matrix based on the solution to the filter algebraic Riccati 
equation, and 

e2) calculating the state estimator gain based on the covariance matrix. 



The process of claim 1, wherein step (e) is performed by 

(el) calculating a covariance matrix E fi"om a solution of the filter algebraic Riccati 
equation in the form 

/ r. \ / \ 

-Z(0 /)0"*(O /)s = o, 

where M is the inertia matrix, 0 is the intemal disturbance matrix and E is the 
external disturbance matrix, and 
(e2) calculating the state estimator gain H from H = S(/ O)'0"* . 
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-14- 



The process of claim 3, wherein step (f) is performed by 



fl) solving H = 



for B and K, and 



f2) setting the optimal mount damping matrix to B and the optimal mount stiffness 
matrix to K. 



The process of claim 1 wherein the state estimator is a Kalman filter. 
The process of claim 5, wherein step e) includes, 

el) calculating a co variance matrix based on the solution to the filter algebraic Riccati 
equation, and 

e2) calculating the Kalman filter gain based on the covariance matrix and the inertia 
matrix. 



The process of claim 6, wherein step (e) is performed by 

(el) calculating a covariance matrix 2 fi*om a solution of the filter algebraic Riccati 
equation in the form 

/ \ / \ 

Z{Q I^@-'{0 7)2 = 0, 

where M is the inertia matrix, 0 is the internal disturbance matrix and S is the 
external disturbance matrix, and 
(e2) calculating the Kabnan filter gain H from = ^{^ O)©"* . 
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8. (original) The process of claim 7, wherein step (f) is performed by 



fl) solving H = 



for B and K, and 



f2) setting the optimal mount damping matrix to B and the optimal mount stiffiiess 
matrix to K. 
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Appendiy D- RelateH Prnceerlingg A ppenHiY 
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Abstract: 

In industrial motor drive systems, a shaft torsional vibration is often generated when a 
motor and a load are connected with a flexible shaft. This paper treats the vibration 
suppression control of such a system. In this paper, a control system design method 
using linear matrix inequality (LMI), which is a tool for control design that replaces or 
complements Lyapunov-Riccati equations, is provided and the speed controller for an 
induction motor by LMI is proposed. In the speed controller, weights are used to 
satisfy tracking and disturbance rejection. Experimental results show the validity of the 
proposed speed controller by LMI, and this controller is compared with the state 
feedback controller 
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Abstract - Li the industiiai motor drive system, a 
shaft torsional vibration is often generated when a 
motor and a load are connected with a flexible shaft. 
This paper treats the vibration suppression control of 
such a system. In thil paiiU', a louUOI systeiri design 
method using Linear Matrix Inequality (LMI), which 
is a tool for control design that replaces or 
complements Lyapunov-Riccati equat ions, is provided 
and the H« spSd tOntrouer lor an induction motor 
by LMI is proposed. In the H<x. speed controller, 
weights are used to satisfy tracking and disturbance 
rejection. Experimental results show the validity of 
the proposed H-> speed controller by LMI. and this 
controller is compared with the state feedback 
controller. 

1. Introduction 

In the industrial motor drive systems such as 
industrial robots, rolling mills, elevators and so on, 
the shaft torsional vibration is often generated when 
a motor and a load are connected with a flexible 
shaft Such a system is called as the two-mass 
resonant system and this paper treats the vibration 
suppression control of this system. 

To suppress the shaft torsional vibration, the 
resonant ratio control, the state fepHhack_control. the ., 
LQG control and the Hoo control are studiedo~6] It 
is often impossible to detect the load speed, the 
shaft torsional torque and the load disturbance 
torque in the real systems. Recently, it is studied 
the new method which estimates unknown state 
variables by using the reduced order observer and 
feed back these state variables by using the pole 
placement design method.[7] But this design method 
doesn't satisfy the fast speed response property and 
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the attenuation property of load disturbance at the 
same time. 

In recent studies, there is the Hoo control which 
satisfies the fast speed response property and the 
attenuation property of load disturbance at the same 
time by using the weients. The H» control is often 
applied in the industrial systems. 

hi recent, because it is possible to represent the 
control specifications and the robust stable conditons 
by Linear Matrix Inequality (LMI), the control 
problem can be transformed into the convex 
problem-fHs] LMIs have gained much attention as a 
tool for control design that replaces or complements 
Lyapunov-Riccati equations. 

hi this paper, the controller design method by 
LMI is provided and the H« speed controller of an 
induction motor by LMI is proposed. Ihe validity of 
the proposed controller is comfinmed by 
experimental results and this controller is compared 
with the conventional state feedback controller using 
the pole placement method,[7i 

2. The two-mass resonant system 

The two-mass resonant system which has a 
motor and a load connected with a fipvihle shaft is 
shown in fig. 1. ihe state equation of the 
two-mass resonant system is as" follows 
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Fig.! The plant model of the two-mass 
resonant system 



JL 



-t 0 



Jl 
0 



(1) 



where, 

Im > Jl ' the motor inertia, the load inertia 

r Bi ■ the motor viscous damping coefficient, 

the load viscous damping coefficient 
KsH '■ the shaft stiffness 
<t>M * 'the motor speed, the load speed 
Tm , TsH , Ti. : the motor torque, the shaft 

torsional torque, the load disturbance torque 

The goals for the control of the two-mass 
resonant system are stated in the following forms, 

(1) Suppress the shaft torsional vibration. 

(2) Inject the effect of the load disturbance torque. 

(3) For the change of the speed reference, the load 

speed follows it quickly without overshoot. 



PLANT P(s) 



-rK>- 

'•(CO-) 



KUCMZKTZD PLANT 



Fig. 2 The generalized plant of the Ht» control 

Tt is the disturbance torque and w is the speed 
command. is weight for rx,and is weight 
for n. If weights are defined as follows, 
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(3) 



then the state space equation for the generalized 
plant and the ratio of output to input, G , are as 
follows. 



(4) 



' X 




\A Bi B-A 




' X 


z 








w 


y 




C2 Z?21 0 




, u , 



A Bi B2 
Cl 0 Z)i2 



At, 0 BrfCw 0 B^Dw Bp 
0 ^ffi, 0 Bttfe 0 0 
0 0 Aw 0 
a 0 0 0 
0 0 0 0 



Bw 0 

0 0 
0 

0 0 



(5) 



3. The H«» generalized plant 

Fig.' 2 is the H« generalized plant to apply the H 
00 control theory to the tv\'o-mass system. In fig. 2, 
toi^ to-i are inputs for disturbance and command, 
respectively, and zi.Zu are ouQDuts. zi is as 
follows. 
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Above G is obtained by choice of W^i , , 
and the Hm controller K is solved by Matlab. 

4. The control system design by LMI 

In this chapter, we represent G of the 
generalized plant by I-MI. Convexity of LMI and 
LMI representation of Riccati inequality for the 
control design are given and convex problems 
involving LMI are treated. 
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Riccati inequality' applying to control system ^ r^V-F'^VF>Q; V: ^ W'^Wy 0 (10) 

where W : scaling matrix, / : level 



design is as follows. 
In Eq. (6), A and B are 



following problems can be obtained. 



Ap 0 B/^^ 

0 Ay^ 0 

0 0 A^^ 



B^[B, 82] = 



0 BJ)^ B, 
Bm 0 0 
0 Bm 0 



(7) 



in real two-mass system and Q , i? , S can be 
sellected by free choice. Eq. (6) can be transformed 
into 



by using Schur complementdoj And Eq. (8) can be 
converted to following general form of LMI. 



(9) 



where 



-A'^Si- SiA SiB] 
B^'Si 0 J 



e, : vector 



Convex problems for solving F{e) > 0 are 
remarkably two problems as follows. 

1) Convex FeasibiHty Problem (CFP) 
find ee??"* such that F(e) > 0 

2) Convex Minimization Problem (CMP) 
minimize c^e subject to F{e) > 0 , 

where is given. 

And because of scaled matrix norm condition, 
II Pi^FH^^II < can be transformed into eq. (10), 

II WFW'W < r <^ yr'^ w'^Vi'F < 7^1 



1) CFP : when matrix F and level r>0 are 
given, problem solving scaling matrix W{ or 

V ): 

^ V: = W'^W>Q such that r^V- F'^VF>Q 

2) CMP : when matrix F and scaling matrix W 
are given, problem minimizing level 7; 

minimize r> 0 such that r^V ~ F'^VF>0 

We can get PF( or V ) in case of minimum y 
by iteration. 



5. H« speed control of the two-mass system 
by LMI 



We design the H« controller of the two-mass 
system by LMI. In fig. 3, the weighting function 
is selected as Pr„ = 10 which does not 
influence the design of controller. The choice of 
weighting functions is very critical because the 
control performances are determined by weighting 
functions. The weighting function is for the 

closed-loop transfer function from the disturbance 
torque Tt to 2i. In order to achieve the design 
specifications such as vibration sujpession and 
disturbance rejection, I I is chosen as 

closed-loop gain plot of the one-mass system. The 
weighting function Wn is for the closed-loop 
transfer function from the disturbance torque n to 
zi , where n corresponds to the speed command or 
noise. The closed-loop transfer function from n to 
q)m is a complementary sensitivity function, so 
W„ is chosen as a high pass filter considering the 
robust stability. Then and W„ are chosen as 
follows. 
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Fig. 3 Bode diagram of weights 
( ver : [dB], hor : [ oj ] ) 



' ^ ' (s+10-')(5+10') ' 

r„ (^^^p5) (11) 

TTie H» controller K{s) can be obtained by using 
Matiab after adjusting the parameters in weighting 
functions as 7^=5, Wrf=40,rr. = 500. 

Weighting functions are shown in fig. 3. The Hen 
controller K{s) is given by 

i^s) = 4.921 X 10^ X 

(5+lxi0M(54-0.6±;80y5+.q.4n 

(5 + 5.6x10^ )(s + 3.59xlO*)(5+9±;33.52)(s-}-0.0011) 

(12) 

6. Experimental system 

The configuration of the experimental system for 
the two-mass resonant system is shown in fig. 4 
There is a flexible shaft ( 1 : 800 [mm], 4> ' 13 
[mm] ) between a motor and a load generator. This 
system is used variable inertia to change the inertia 
of a motor or a load generator. For this paper, it is 
only used variable inertia for a load generator. The 
mechanical constants of the two-mass resonant 
system are listed on Table 1. 




Fig. 4 The experimental system for the two-nnass 
resonant system 



Table 1 The nr)echanical constants of the two-mass 
resonant system 



Jm 


Jl 


^Sff 


0.008 [kgm'] 


0.08 [kgm'] 


50.527 [Nm/rad] 



7. Experimental results 

Experimental results show the validity of the 
proposed controller and this controller is compared 
with the state feedback controller. The inertia ratio 
is 0.008 : 0.08 ( : Ji ) in tiiis paper 

Experiments are composed of speed response 
experiment and load disturbance response 
experiment. 

7.1. Speed response experiment 

The reference speed of the load is 200[rpm] as a 
step input at 0 [sec]. Fig, 5, 6 are experimental 
results of step response under the state feedback 
control. Fig. 5 shows load speed and motor speed 
and fig. 6 shows shaft torsional torque. The speed 
transient time from 0 to 200 [rpm] is about 0.7 
[sec]. If the transient time becomes shorter, 
overshoot of the load speed is Iai:ger. In the 
transient time, maximiim tongue is 6 [Nm]. 
Fig. 7, 8 are experimental results for control by 
UVn. Because the speed transient time from 0 to 
200 [nml is about 05 [sec], the transient tinrie of 
proposed LMI controller has much less than that of 
state feedback controller. This is because of 
adjusting the weighting function for speed response. 
In the transient time, maximum torque is 9 [Nm]. 
But this maximum torque does not influence the 
two-mass system 
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Fig. 5 State feedback controller 
150 [rpm]/div, 0.5 sec/div 
(upper) Real motor speed, (down) Real load speed 



Fig. 9 State feedback controller 
150 [rpmVdiv. 0.5 sec/div 
(upper) Real motor speed, (down) Real load speed 





Fig. 6 State feedback controller 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 



Fig. 10 State feedback controller 
3 [Nm]/div, 0.5 sec/div 
Real shaft torsional torque 





Fig. 7 LMI controller Fig. 11 LMI controller 

150 [rpml/div, 0.5 sec/div 150 [rpm]/div, 0.5 sec/div 

(upper) Real motor speed, (down) Real load speed (upper) Real motor speed, (down) Real load speed 





Fig. 8 LMI controller 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 



Fig. 12 . LMI controUer 
3 [Nml/div, 0.5 sec/div 
Real shaft torsional torque 
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7.2. Lxtad disturbance response experiment 

A resistor load which can be connected with DC 
generator is used as a load. By using a resistor 
load, the load disturbancs torque is increased from 0 
to 4.5 [Nm] at 200 [rpm]. 45 [Nm] is 37% of rated 
torque. Fig. 9, 10 are experimental results for the 
state feedback control. Load speed dropped 17% of 
200 [rpm] and recovered after 1.2 [sec]. This is 
because the inertia ratio is big. If the inertia ratio 
increases, the recovery time is much longer. 

Fig. 11, 12 are experimental results for control by 
LMI. The load speed dropped 5% of 200 [rpm] and 
recovered after 0.6 [sec]. Then, the speed drop of 
the proposed LMI controller has much less than that 
of state feedback controller. The recovery time of 
the proposed LMI controller has much less than that 
of the state feedback controller. 



8. Conclusion 

This paper treated the speed control of the 
induction motor system which has a motor and a 
load connected with a flexible shaft The H« speed 
control of an induction motor with two-mass 
resonant system by LMI is proposed. In the 
proposed control, the H« control, which satisfies the 
fast speed response property and the attenuation 
property of load disturbance at the same time by 
using the weights, is used and is r^Desented by 
LMI which is a tool for control design that replaces 
or complements Lyapunov-Riccati equations. 

The experiments show the excellent performance 
of the proposed Hoo speed control by LMI compared 
with the state feedback control. 
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"Optimum Vibration Control of Distributed Parameter Rotor Bearing System by Using 
j Electromagnetic Bearing" 
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Abstract 



[Associate Professor, Hankuk Aviation University 

As the modern rotor systems have a tendency of being operated at high speed and 
designed to be lighter, the control of excessive vibration becomes more and more 
important because of increasing requirement In accuracy and reliability. Since the 
flexible rotor bearing system is a typical hyperbolic type distributed parameter 
system, the control, passive or active, not accounting for such a nature often does 
not yield satisfactory performance. The primary objective of this study is to develop 
[efficient control design methods which account for and use the distributed parameter 
and isotropic properties. The output feedback control scheme accounts for the 
practical aspects such as the observation and related stability problems caused by 
spillover effects, while the optimal design deals mainly with the theoretical and 
numerical aspects. Optimal pole assignment in distributed parameter control systems 
I is established and the existeRee-of optinqal solution is proven in case of compact 
j linear feedback. The^ntethod, an extensidn of the finite dimensional design, gives the 
I solution of the operator Riccati equation/which is known to be one of the most 
j difficult problems fa6^by the designpf^. By adopting the complex notation for the 
isotropic systems, it is^howuxlhat-tHe closed loop eigen structure is determined, 
contrary to the real field design, so as to preserve the isotropic property. In order to 
utilize the xymmetric property between the two directional vibrations: vertical and 
horizontal or forward and backward, a new control scheme is designed by introducing 
the complex state for the systems having the isotropic property. The algebraic 
relations concerning distribution of the eigenvalues of complex matrices are 
developed so as to duster the poles into specified regions, rather than fixed positions, 
imposing the relative stability margin, e.g.. uniform damping or uniform damping 
ratio. The closed loop eigenstructure has the direct one-to-one relation with the 
physical parameters and the isotropic property is preserved under feedback control. 
[The method transforms Riccati equation to Lyapunov equation of the same dimension 
so that computational effort is almost removed ^ It is shown tb atthc method can be 
^ used effectively in vibration control of rotating mac hinery such as the isotropic rotor 
tbcaiUiU ^ystem, i:"he magnetic bearin g system an JTT re lu l a L i ii g LHCU l j r-4igkJo 
l^achte^e thg ^uQl objuLliv c o f i ^l dUill/a tion and spillovei buppiei>bT Oftre-^f>strained 
output feedback control scheme is developed based on the optimization of a modified 
performance index which includes spillover terms. In order to construct the 
appropriate performance index, the method makes use of the two-time scale property 
of the primary and secondary modes. The control model is constructed. based on the 
singulariy perturbed modal model, which retains the simple modal structure while 
preserving the accuracy of the frequency domain modeling. A set of optlmality 
condition is derived and iterative solution procedure is also, presented. A set of 
simulation examples show the effectiveness of the proposed scheme. In order to 
control the rotation related periodic disturbances, a disturbance acconlmodating 
controller is also designed based on the disturbance estimator. The. main, ad vantage 
I of the proposed scheme is the simplicity in numerical calculation and estimator 
structures, requiring only a solution of non-symmetric Lyapunov equation of relatively 
small order. Lyapunov equation is the first order perturbation solution of the matrix 
! cubic equation which yields the exact solution. Through the rnatrix perturbation 
analysis, it is proven that, for small gain, the approximate soUitioh converges to the 
exact one. The control scheme is proven to effectively, control the build up vibration 
I subject to sudden imbalance through the numerical siiftUiatiohai. Finally, the output 
feedback control experiment is performed by using a magnetic bearing when the rotor 
system is subject to external impulse and imbalance. The magnetic bearing with four 
[magnets is manufactured and the analog design is transformed to.the.discretized 
form and implemented by using a 16 bit microcomputer. The.contrbl model. is 
constructed by using the experimentally obtained modal parameters through a set of 
modal testings in complex coordinates. The proposedxontror'scheme .enables the 
rotor to run through the first.critical speed without excessive Vibration against the 
Initial system unbalance. In order to verify the performance of the proposed control 
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initial system unbalance. In order to verify the performance of the proposed control 
scheme, the. results are compared with the numerical ones. It is shown that the 
experimental results are fairiy well coincident with the predicted ones. In conclusion, 
through the extensive analysis and experiments, it is shown that the proposed 
controllers can be used effectively in practical implementation as well as in theoretical 
and numerical aspect of controller design for the rotor bearing systems. 

"I 2004 Vibration Control Lab. All rights reserved. 
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Abstract: " ' 

The design, development, and application of a research servowriter, which is necessary 
to record magnetic position encoding schemes, are presented. The ^ervowriteMs ^ 
J ?he"d?.. Jh?!;"'^ ^""I'^'^^l djsk drive technologies. Submicron positioning accuracy 
o the disk actuator is achieved by integrating the state-of-the-art electronic, optic 
electromechanical, and digital signal processing hardware. The servowriter alLo provides 

~or J^^^^^^ f ^^^^-^^ ^ programmable pattern " 

generator and has several modules of customized electronic hardware, each controlled hv 

lsl7r:ZTrrZTT. '^'^ ''^^ (read/write) head and actuator ' 

assembly IS measi^red by a laser mterferometric system. The measured position is used 

vaSbTtUd ^^"^ ^^'^'^ "^-d to implemeSra state- ' 

variable-based proport.onal-mtegral-derivative controller on the Intel (8086/8087) 
microprocessor chip set. Both systematic and random error sources associated with 
seifm\'l:Sand f '^"f ^nvironmentai™!^^^^^^^ 

^ oS arf idenNJ^^d .^H^'f '^l"''" ^^'^'^ '"^""^^^'y ^°ntnbute to servowriting 

errors, are identified, and steps to contain these error sources are described 
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Development of a Servowriter for Magnetic Disk 

Storage Applications 

DANA H. BROWN and MUTHUTHAMBY SRI-JAYANTHA, member, lEEr. 



Abstract— Tht design, development, and application of a research 
servowriter^ which is necessary to record magnetic position encoding 
schemes, are presented. The servowriter is designed to evaluate emerg- 
ing disk drive technologies. Suhmlcron positioning accuracy of the disk 
actuator is achieved by integrating the state-of-the-art electronic^ op- 
tic, electromechanical, and digital signal processing hardware. The 
servowriter also provides flexible and sophisticated signal encoding 
schemes through a programmable * 'pattern generator'* and has sev- 
eral modules of customized electronic hardware, each controlled by a 
supervisory computer. The position oflbe R/W (read/write) head and 
actuator assembly is measured by a laser inlvrferomelric system. The 
racusured position is used to estimate the velocity of the R/W head, 
which is then used to implement a stutc-vartable-based proportional- 
integrat-derivative controller on the Intel (8086/8087) microprocessor 
chip set. Both systematic and random error sources associated with 
phase and amplitude encoding schemes are described. Environmental 
factors, such as seismic vibration and temperature change, which in- 
directly contribute to servowriting errors, are identified and steps to 
contain these error sources are described. 



1, Introduction 

THE HIGH -performance disk drives of loday arc pre- 
cision rotating machines, with head positioning 
mechanisms capable of achieving submicron track fol- 
lowing accuracy. As the demand for track density in- 
creases, the ability to provide a higher resolution position 
signal will become mandatory. In order to maintain the 
magnetic R/W (read/write) head position along the center 
line of a desired track, a magnetic pattern that indicates 
the position of the R/W head is prerecorded on the disk 
surfaces during the manufacture of the disk drive. The 
device used to record the magnetic pattern is called a ser- 
vowriter n, P- 719]. Since the magnetic pattern helps the 
servo control algorithms perform the track-following or 
track-accessing task, it is also referred to as a ser\'o pat- 
tern. 

Of the several major steps involved in the manufacture 
of a disk file, the servowriting operation takes place just 
prior to sealing the head/disk assembly in a clean room 
environment. A servowriter in itself is a complex device 
that integrates state-of-the-art electronic, optic, and elec- 
tromechanical hardware with sophisticated signal encod- 
ing schemes to achieve the submicron accuracies that are 

Manuscript received March 6, 198**; revised November 6. 1989. 

D. H. Brown is with IBM Applications Business Systems. Dcpanmcni 
215, Rochester. MN 55901. 

M. Sri-Jayantha is with the IBM Thomas J. Watson Research Division. 
Yorkiown Heights, NY 10598. 

IHHE Log Number 8933409. 



demanded in the industr>'. There are three broad catego- 
ries of position encoding techniques used by the disk drive 
industry: 

1) nonmagnetic position encoding (optical sensors 
mounted on the head actuator, mechanical detents 
for each track, stepper motors, etc); 
.2) magnetic servo pattern based position encoding 
(amplitude encoded schemes [2] and phase encoded 
schemes [3]); 

3) combinations of the above (optical sensors with a 
single servo burst per track to remove low frequency 
track- following errors). 

The servo control of the IBM9332, for instance, is based 
on a phase encoded servo pattern [4], whereas the 
fBM33XX series is based on an amplitude encoded servo 
pattern. A short overview of servo patterns is given in the 
Appendix. 

This paper addresses the design, development, and ap- 
plication of a research servowriter that is targeted to 
achieve the magnetically based position encoding tech- 
nique in category (2). Following an overview of the ser- 
vowriter, its major elements are described briefly. The 
servowriter electronics are made up of several customized 
electronic hardware modules, including an Intel (8086/ 
8087) microprocessor chip set. Each module or subsys- 
tem is associated with a function that is required to sup- 
port the ser\'o writing and testing process and is controlled 
by a supervisory computer (JBM-PC/AT), The pattern 
generator is software configurable and is controlled by the 
IBM-PC/ AT. The JBM-PC/AT is menu driven and it is 
through this menu that the user controls or monitors the 
servowriting operation. More details are presented in |5|. 

The position of the R/W head is measured bv a laser 
interferometric system [6], The laser based position is 
used to estimate the state variables of the. R/W head oo- 
silioning mechanism. These estimates are then used to . 
implement a slate-variable-based proponional-integral- 
derivative (PID) controller, which positions the mecha- 
nism with zero steady -state error in the presence of bias 
force. Estimator/controller algorithms are implemented in 
8087-based floating-point arithmetic using assembly lan- 
guage. An example of a closed-loop transfer function of 
a servowriting mechanism is presented, and factors that 
would contribute to track-following error, which is det- 
rimental to the servowriter performance, are discussed. 
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Details related to clock pulse generation that is man- 
datory for locating the pattern generator pulses are de- 
scribed. The servo pattern is generated from 1024 X 32 
bit memory-based data, thus providing flexibility to study 
the performance of various encoding schemes. Finally, 
the error sources associated with encoding schemes, and 
environmental factors <^iich as seismic vibration and hur 
"^'^itv. which indirectiv contribute m <:''rvowriiiBe crrorx. 
are identified and steps to contain these error sonrcc^i are 
aescnbeo. 

n. Servowriter DcSCRfPTrON 
The servo patterns are written on one or more disk sur- 
faces, depending on the design of the file being written 
(i.e., dedicated versus sector servo). The servo patterns 
can be written with heads mounted on the disk drive ac- 
tuator, or a special set of heads can be mounted to the 
servowriier to perform this task. The servowriting process 
is a repetitive process of writing a track and moving the 
head to the next track. The process repeats itself until all 
tracks are written. The process of writing the servo pat- 
terns is relatively straightforward, but since it must be 
precisely done, the circuitry to do this can become quite 
^. complex. Fig. J shows the major elements of a servowri- 
tcr, some of which are briefly described now. 

The heads must have a sensor mounted to them to in- 
dicate their radial position on the disk. Typically a laser 
interferometer is used to measure the head position. This 
sensor system is used in a closed-loop servo system to 
locate each track to be servowritten. Once the head is po- 
sitioned over a track, the pattern is written. However, the 
pattern must be synchronized to the disk rotational posi- 
tion so that adjacent concentric tracks line up according 
to a selected encoding rule. For this task, a sensor is used 
to supply fine resolution clock pulses that are synchro- 
nixed to the angular position of tlie disk. This sensor is 
typically a stationary RAV head called a clock head. As 
the spindle motor speed varies due to random distur- 
bances, the clock pulses follow the speed changes accord- 
ingly. 

An encoder or pattern generator supplies the correct 
servo code to the heads for writing on the disk. This pat- 
tern generator follows the angular position clock as it 
writes each track so that the servo pattern on that track is 
aligned to the servo pattern that was written on the pre- 
vious adjacent track. The pattern- generator is pro- 
grammed to change the code written on each track so that 
the servo signal on read back will indicate over which 
track the head is positioned. 

Each subsystem, together with the tasks performed by 
It, is coordinated by a central supervisory computer. A 
suitable architecture is used to share the computing, con- 
trol, and monitoring tasks required by a servowriter. 

III. Major Elements of a Servowriter 
A, Computer Archiieciure 

Fig. 2 shows the microprocessor hardware interconnec- 
tion necessary to achieve the servowriting operation. The 
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Fig. 2. Computer architecture. 

servowriter is controlled by an IBM PC/AT with external 
data output (DO) and data input (DI) registers. The con- 
trol program for this servowriter was written in IBM in- 
terpreter BASIC with some data manipulation programs 
(such as the pattern generator down-loader) written in 
8088 Assembly language. The memory in the PC/AT is 
640 kilobytes. However, the program to control the ser- 
vowriter is conuined in a 64 kilobyte segment. Ahhough 
interpreter BASIC is slow compared to Assembly lan- 
guage, 2000 servo half tracks can be written on a disk in 
well under I h. This is considered adequate for a research 
servowriter, but is not acceptable for a production ser- 
vowriter. The software would have to be optimized to re- 
duce that time by a factor of at least 5. 

B, Digital Position Control System 

The head is positioned by mounting a retroreflector to 
the actuator and tracking it using laser interferometry, A 
closed-loop digital servo feedback system is employed to 
hold the head position to vC'ithin a few microinches of the 
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target position. A second-order actuator model with an 
augmented position integral state is used for contrnller de- 
sign. A moving coil actuator with a current-mode driver 
is thus modeled as follows: 
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where 

c<f damping constant (N • s/m), 

kf? actuator force constant ( N/A ), 

ks stiffness constant ( N / m ) , 

k,. sensor gain constant (laser counts/m), 

m equivalent actuator mass (kgr), 

u control current commanded by servo writer pro- 
cessor (A), 

V laser measurement noise (equivalent laser count), 

process noise vector (appropriate units), 

X ^ f position ( m ) . integral ( m - s ) , velocity ( m /s ) ] , 

.CRcr reference position (m), 

I measured position (laser counts). 

The representative model defined by (I) is used to de- 
sign a PID-type controller using a pole-placement design 
technique. The introduction of process and measurement 
noise terms emphasizes that the system is not determin- 
istic, and that the tracking performance, therefore, can 
never be ideal (i.e., zero error). The unmodeled distur- 
bances act through the orocess noise term, H!. 3nd piea- 
surement noi.<ie annears through the noise term, v, seis- 
mic tor ground) vioraiion not only airects tne accuracy of 
uic inouei aennea oy ( i ). but also contnt>utes to measure^ 
meni noise bv causinn the laser optics to move relative to 
ui6 actuator. Isolation ot these sources contributes signif- 
icantly to the cost of a servowritcr. 

The continuous model is discretized for controller/es- 
timator design. A sampling rate of 300 /zs is used to dis- 
cretize the model and is represented as follows (ignoring 
noise terms, and setting x^cr = 0): 

x^., = 4>x, + Tm, (3) 
= h^x, (4) 

where 
k sample index, 

h measurement matrix as defined in (2). 

Equations (3) and (4) are used to select the controller 
feedback gain constants. A third-order Butterworth pole 
pattern is used to locate the desired closed-loop poles of 
the position servo. 



Since the measured position (i.e., laser count) is not a 
continuous variable, a reduced order "current" estimator 
is designed to obtain the velocity estimate [7, p. 149). In 
the design of the estimator, measured position is assumed 
to be exact, and the position (error) integral state is ex- 
cluded from the estimator model. Therefore, the velocity 
estimator is of first order, and the pole-placement design 
yields a single gain corresponding to a desired estimator 
real pole location. The integral state is derived using the 
trapezoidal rule for integration. 

The controller/estimator gains and the system matrix 
elements are arranged together so that each laser position 
measurement is used to generate an appropriate electric 
current command in an explicit state variable fonn. The 
controller/estimator computations are realized in floating 
point arithmetic. Using a Hewlett Packard 3562A dy- 
namic signal analyzer, a closed-loop transfer function 
(CLTF) of an experimental sysieni was obtained. Fig. 3 
shows the measured CLTF of an integrated disk/actuator 
system under servowritcr control. The bandwidth of this 
system is limited by a 400-Hz resonance mode specific to 
the actuator used (in which the mounting of the retrore- 
flector accentuated the 400-Hz mode). The ser\'owritcr 
controller/estimator algorithm has adequate computa- 
tional speed to support a closed-loop system with much 
higher bandwidth if the actuator modes are not a problem 
from a stability standpoint. The presence of the integral 
state provides the required low-frequency disturbance re- 
jection. In order to suppress high-frequency disturbance, 
greater bandwidth with minimum CLTF peaking is rec- 
ommended. 

C. System liming and Clock Writer Design 

There are different methods in practice today for writ- 
ing clock signals. The simplest of all methods is to let the 
clock writer run freely at a designed clock frequency and 
to dither the disk speed (because It can never be regulated 
perfectly at a constant speed) until closure (i.e., an inte- 
gral number of complete clock cycles around a track) is 
achieved. This task is automatically performed by the 
clock writer in 5 to 30 s. Hence, as defined in the next 
paragraph, the total number of clock cycles is first se- 
lected for a given problem, and through repealed write 
trials, the clock signal is recorded until satisfactory clo- 
sure is met. 

The total number of clock cycles per track is based on 
four parameters that have various system implications. 
They are defined as follows: 

A servo pattern sampling period (s), 

^'rpm spindle motor speed (rpm), 

/clock clock track frequency (MHz), 

"Section integer number of sectors per revolution. 

The sampling time is specified to meet the performance 
requirements of the disk actuator servo. (Note that the ser- 
vowritcr laser sampling rate is independent of the disk 
servo sampling rate.) The spindle speed is usually deter- 
mined from magnetics considerations and is set to a spec- 
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Fig. 3. Closed loop irdnsl'er funciion of an experimental system. 

ified value. The clock frequency should meet the R/W 
head frequency characteristics as well as the requirements 
imposed by ihc specific servo encoding scheme. The 
number of sectors is. however, constrained by the rela- 
tionship 
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The initial estimate of /w.ccrnr will not necessarily be an 
integer. Therefore, a convenient integral number of sec- 
tors, /w*aorr corresponding to an exact, but slightly ad- 
justed, sampling rate A**, is selected as the first step. For 
instance, in order to drive the brushless dc spindle motor, 
an even number of sectors divisible by 12 (related to the 
number of motor poles) may be necessary. Based on the 
design requirements of the servo pattern, n^,^a, clock 
cycles at a frequency /,,,^,^ are to be generated bet^fveen 
two consecutive sectors with a duration dictated by an ad- 
justed sampling rate A*. This requirement is met by the 
relationship 



(6) 



In the case where Hsccior of (6) does not result in an integer, 
it may be rounded to the nearest integral value » /i^^^tor 
Fractional adjustment of clock frequency (in megahertz) 
would typically satisfy this constraint. Hence, the total 
number of clock cycles per revolution of the disk spindle 
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This redefines the spindle speed to be 
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(7) 



(8) 



which, though not exact, is near the desired spindle speed. 
For example, a consistent set of parameters would be A* 
= 99.206 /xs, = 248 [=99.206 X 10'^ x (2.5 X 
10*)]. "toi.! = 248 X 168 = 41 664 clock cycles per rev- 
olution, and ^;^p = 3600.23 ipm, where the clock track 
writer is programmed to write n^o^, (=41 664) clock 
cycles at a frequency of 2.5 MHz. This is done by turning 
on the write current to the clock head while writing 2.5- 
MHz data to it. In order to have exactly 41 664 clocks 
written under the clock head, a counter must register how 



many cycles have been written during each revolution of 
the disk. When exactly 41 664 clocks have been written 
for one revolution (achieved by dithering the speed of the 
motor slightly around A^rpm). the write current is shut off. 

Counting clock pulses during clock writing requires a 
reference mark to indicate precisely when the disk has 
turned through a single revolution. This is done by writ- 
ing an index mark with the servo head before clock writ- 
ing. The index decoder (see Fig. 1) monitors the read- 
back signal from the servo head and outputs a pulse each 
lime die index mark is detected. This pulse (called index) 
initially starts the clock writer counter during clock writ- 
ing and resets it after each unsuccessful pass of clock 
writing. Clock track verification is done by simply read- 
ing back the written clock track and counting the number 
of clocks in one revolution. If disk defects damage one or 
more clock pulses, this will be reflected in the final count. 
An index mark is also encoded in the clock track to in- 
dicate the start of the clock track. Three cycles of 1.25- 
MHz frequency marker can be embedded in the clock 
track without seriously affecting the 2.5-MHz clock sig- 
nal. This index is used to synchronize the pattern gener- 
ators at the start of each servo track. 

The clock track resolution of 400 ns is too coarse for 
writing servo patterns. A type 2 phase-locked loop (PLL) 
that multiplies the clock track frequency by 16 and holds 
the output phase of this 40-MHz signal constant with re- 
spect to the clock track is employed. The PLL output ef- 
fectively subdivides the disk into 666 624 equally spaced 
lime slots of 25 ns each. 

D. Motor Speed Conirol 

The speed control of the disk is critical during clock 
writing and clock read-back. If it varies too much during 
either clock writing or read-back, good closure will not 
be attained and the PLL that generates the system time 
base will go out of lock. This will make the writing of 
good servo patterns impossible. The speed of the disk is 
controlled by comparing motor Hall sensors (or, more ac- 
curately, the clock track signal) to a stable reference 
clock. For a research scrvowriter, the reference is derived 
from a high-precision (six to seven significant digits) fre- 
quency source. This is necessary for changing the number 
of clock cycles by a couple of counts per revolution. For 
5.25- or 3.5-in disk drives, the motors are small enough 
in power consumption that single IC motor controllers are 
readily available, and with the proper compensation, they 
can control the motor speed to about 0.01 percent of ref- 
erence speed. 

E. Pattern Generator 

The circuitry that writes the servo patterns on the tracks 
is called a pattern generator. It is usually programmable, 
but for servowriters with simple servo patterns it can be 
hardwired to write a limited number of servo codes. The 
pattern generator is a clock track counter that begins 
counting the 40-MH2 PLL output clocks when the clock 
track index is received. At certain preprogrammed counts, 
the output of the pattern generator changes states to re- 
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verse the write current in the servo head. It is this se- 
quence of transitions that makes up the scr\'o code. 

There are pattern generators of differing designs on the 
commercial test equipment market today. The pattern 
generator designed for a typical IBM servowriter is sim- 
ilar to some of the commercial ones except for instruction 
modifications to facilitate writing servo patterns. The out- 
put of the pattern generator is a shift register that clocks 
out the data with each 25-ns clock. The shift register can 
be loaded with 9-16 bits with each pattern generator in- 
struction and then shifted out at the system PLL clock rate 
(40 MHz). 

The pattern generator contains a block of memory, for 
example, 1024 x 32 bits. Each 32 bit word contains 3 
bits of instruction op code, 3 bits to indicate the shift reg- 
ister data length, 16 bits of data, and 10 bits of op code 
modifier (i.e., subroutine call address, branch addresses, 
or loop count information). Although the in.structions are 
simple, they allow the writing of any servo pattern used 
to date and, in most cases, all tracks for an entire disk can 
easily be contained in the 1024 locations. 

F. Pattern Testing 

Writing servo patterns on a disk at exactly the correct 
locations does not guarantee a correctly written disk. It is 
possible that defects on the surface of the disk were pre- 
sent and that the servo information was not recorded at 
those points. The recording head may have had a shorted 
turn or other defect that prevented the servo pattern from 
being written on the disk also. 

The pattern test circuitry allows a read-back check to 
verify that a good servo pattern has been written. While 
the concepts and circuits for pattern testing are not diffi- 
cult, this part of the servowriter is troublesome. Errors 
are neither black nor white. They fall into a broad range 
of values, many landing on the threshold that the pattern 
test uses to evaluate the correctness of the pattern. The 
result is that a particular disk defect or servo error may 
pass on some tests and fail on others, or it may pass on a 
particular servowriter and fail on another. The pattern 
testing spans several subsystems and the details are pre- 
sented in [5]- 

IV. Error Sources 

The servowriter errors can be either random or sy.stcm- 
atic and are the subject of the next section. 

Error Mechanisms in the Servo Pattern 

The phase encoded servo pattern looks at the average 
phase of many transitions with respect to a reference fre- 
quency (crystal). As a result, it is relatively immune to 
single transition defects and small amplitude variations, 
[deally, the servo pattern phases should change continu- 
ously as the head position varies. Unfortunately the head 
gap geometry requires a discrete approximation. This ap- 
proximation to a sloping phase angle is a 5/2 ns discrete 
time shift for every 1/2 track width of position change. 
Because of thh discrete implementation, there can be rip- 



ple in the position signal generated by the demodulator as 
the head is displaced. Selecting a finer discrete resolution 
does not solve the problem of ripple completely. Mag- 
netic heads write tracks that are slightly wider than the 
width of the head gap. Each pass of the head while ser- 
vowriting leaves a noise fringe field between tracks. This 
fringe field is a fixed width that is a bigger percentage of 
each track width as track spacing dimensions shrink. If 
quarter tracks are written with 5/4 ns phase shifts, the 
linearity should get better (i.e., the ripple amplitude 
should decrease), but a larger percentage of the servo pat- 
tern would be made up of fringe field. Thus the improved 
linearity is more than negated by the noise induced. 

In amplitude servo schpmes, the demodulator output 
versus the actual head position becomes nonlinear as the 
head moves away from the track center line. This nonlin- 
earity is very much dependent on the ratio of the head 
width to the track width. Thus the sen'o gain of the de- 
modulator output (millivolts/micrometers) varies signifi- 
cantly with head width. Some amplitude servo patterns 
use only a single transition in each position sample. If 
such a transition is damaged, though, the position signal 
for that servo sample would be incorrect. Longer burst 
servo patterns overcome this problem to a great extent by 
averaging several redundant transitions (as shown in Fig. 
6(b) in the Appendix). However, this redundancy cannot 
be carried too far because longer servo patterns begin to 
cause errors during high-speed seeks. They also consume 
the useful data space in files with sector ser\'o architec- 
tures. 

1) Tangential Delay: Que to the skew of the head on 
the rotary actuator as shQwn in Fig. 4, the tops of the 
transitions on track A' -f 1 do not line up with the bottoms 
of the transitions on track A' as shown in Fig. 4(a). How- 
ever, if the pattern written on track N is delayed slightly, 
they will line up as shown in Fig. 4(b). For such small 
adjustments, the 3-ns resolution is necessary. Besides 
head skew, nonradial motion of the actuator introduces a 
timing shift in the K, or tar^gential, direction, as shown in 
Fig. 5. The tangential delay circuit can be programmed 
to a new correction factor for each radial servo track to 
remove these systematic errors. There is .still a random 
delay due to the uncontrolled vibration of the clock head 
(S^t-uick) and the servo head suspension (SVsen-o)- These 
two random components are time dependent, and are 
shawn alone the .Vand 7 direction for simplicity, tne ef- 
fect ot random vibratjon can be reduced by special treat- 
ment of the head suspensions or by selecting low jitter 
heads. Fig. 5 shows a generalized delay generating mech- 
anism relating the track spacing (6r), radius (r), disk ro- 
tational velocity (cj), and slope of the actuator path 
(dy/dx). The relative lime delay of a second track writ- 
ten next to an ideal first track (i.e., no head vibration) can 
be shown to be 

Relative Delay 

^{SVscrvo - <5.T,.iock) + ~ 5rj 
- [t X 5r] ^ — r- (9) 
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Fig. 4. Tnin.sitions with and without langctiai delay correction. (The edge 
of a transition can either be orthogonal or inclined (as shown) to the track 
ccnlerline.) (a) Without langelial delay correction the transitions do not 
line up track to track, (b) With tangential delay correction the transitions 
do tine up with adjacent tracks. 




Fig. 5. Effect of nonradial head motion on tangential delay. 



where t is the ideal delay per unit distance of r. There- 
fore, the first group of terms in the right-hand side of (9) 
represents the desired ideal time delay. However, in prac- 
tice the parameters defined in the second group of terms 
also contribute to the total delay. Observe that the delay 
parameters become more complicated when the actuator 
travels at higher velocity, as in the case of a seek. 

2) Spindle and Head Suspension Jitter: As pointed out 
above, track to track timing alignment is essential. Me- 
chanical movement of the spindle center line as it rotates 
will cause the disks to shift the track center line contin- 
ually under the servowriter heads, even if the access sys- 
tem holds the heads perfectly stationary. Radial nonre- 
peatable spindle motion can cause the servo pattern being 
written to squeeze toward one of the adjacent tracks. Tan- 
gential nonrepeatable spindle motion can cause the data 
to be written out of alignment with the adjacent track, 
even if the timing signals were to occur at the proper 
times. Head suspension jitter can cause symptoms similar 
to spindle jitter. The head suspension jitter comes from 
both the clock head and the servowriting head. The jitter 
primarily occurs at right angles to the head suspension's 



main axis, as the stiffness of the suspension is weakest in 
that direction. In-line heads reduce tangetial jitter. 

3) Vibration Isolation and Mechanical Alignment: In 
order to minimize file vibration during servowriting, the 
file is isolated from external vibrations by mounting it on 
a large granite base that floats on air isolation mounts. 
Ideally, the file should be clamped at three points, so that 
the base casting is not distorted while servowriting. in 
practice, bending modes of the casting sometimes can oc- 
cur with just three clamps holding the file down and these 
vibrations can be written imo the servo pattern in the form 
of track positioning errors. Thus some experimentation 
will usually have to be done to determine the bcsl clamp- 
ing arrangement for a particular disk drive. 

The alignment of the disk and head stack for sector ser- 
vowritten files is critical. During file operation, all heads 
on an actuator should be centered over the same track 
number on each surface. In addition » the servo patterns 
on all surfaces should pass under the heads at the same 
time. At the lime of servowriting, this alignment is veri- 
fied. However, if the file is clamped incorrectly in the 
servowriter or in the final product configuration, the ac- 
tuator and disk spindle will tilt, resulting in radial mis- 
alignment and timing shifts from head to head. Sector 
servo files are designed to allow for some radial and tim- 
ing misalignment, but typically these numbers are on the 
order of 5 /im and 500 ns. respectively. 

4) Environment: The key to successful servowriting is 
cleanliness and the precise placement of servo infonna- 
tion. For this reason, the environment of the servowriter 
is extremely important. A servowriter should be in a clean 
room so that exposed files will not become contaminated 
with dirt particles. Even in a clean room, though, dirt will 
accumulate on any surface (although minimal). A simple 
cover should be placed over each open file in the clean 
room while it is waiting for an assembly operation or ser- 
vowriting. During the servowriting operation the cover 
can be removed temporarily. Besides keeping dirt off of 
the disks and out of the head-disk a.ssembly (HDA) en- 
closure before it is sealed, magnetic particles should be 
purged from the clean room. Temperature is critical in 
precision placement of the tracks. The file dimensions will 
shift as the temperature changes, so it is best to keep the 
temperature between 68 and 72 °F to minimize these me- 
chanical changes and to facilitate operator comfort. Hu- 
midity should be between 40 and 60 percent. Above 60 
percent, the humidity causes the servowriter to rust. Be- 
low 40 percent, static electricity begins to cause prob- 
lems. 

V. Conclusions 
Servowritcrs may seem complex on the surface because 
many interleaved subsystems are involved. However, each 
subsystem is easy to understand. Becau.se of the building 
block approach, the servowriter is very logically con- 
structed. It is a precision machine that should be treated 
with care. Time and money spent on the development of 
a servowriter usually are recovered in a less expensive 
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product. However, a servowriter cannot be designed and 
developed independently of the file. The servowriter' s 
mechanical clamping and access systems are directly tied 
to the mechanical design of the file. Likewise, a file may 
behave very poorly after servowriting if its design has not 
been thoroughly verified on the servowriter However, no 
matter how much planning and testing are performed, in- 
evitable adjustments to the file and servowriter do occur 
during development. They usually come in the areas of 
vibration isolation, file alignment, and actuator position- 
ing. For this reason, mechanical design considerations in 
the file and servowriter should be given high priority. En- 
hancing either the software or the electronics alone will 
not compensate for mechanical design deficiencies. The 
environment causes problems that are not directly trace- 
able. Magnetic particle contamination may appear once 
the product has left the servowriter. Errors in servowrit- 
ing induced by temperature shifts or head crashes due to 
accumulation of debris never directly point to the cause, 
thus requiring anxious search for solutions, usually later 
during high-volume production. With care, it can be an 
extremely powerful analysis tool, and without it, a never- 
ending source of problems. 

Appendix: Servo Patterns 
The phase encoded servo pattern (or phase modulated 
servo) works on the principle of varying the relative phase 
of two bursts of squares waves (both at the same fre- 
quency) for each successive track. As a head crosses 
tracks, the phase shift between burst 1 and burst 2 for each 
track changes. This phase shift becomes an indicator of 
track position. In Fig. 6(a) the vertical dashed line rep- 
resents the boundary between burst 1 and burst 2. When 
the head is positioned over track + I, there will be no 
phase shift between burst I and burst 2. When the head 
moves to track N 2, there is a 90° phase shift between 
burst I and burst 2. The phase shift varies quite linearly 
with position for all head locations between tracks + I 
and A/ -f- 2. Note that the pattern repeats every four tracks 
so that the same phase shift .results when the head is over 
track A^, + 4, + 8, • • • . The benefits of a phase 
encoded pattern are that the servo signal frequency does 
not vary, the pattern is linear over four tracks, and single 
bit disk defects usually do not have a significant effect on 
the servo position signal, since multiple cycles are aver- 
aged. 

The amplitude servo scheme is more commonly em- 
ployed today. It has good linearity over track centers, but 
becomes nonlinear as the head moves off track center by 
more than 35 percent of the track width. This scheme is 
usually easier to demodulate and to servowrite. Ampli- 
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tude servo relics on the principle that signal amplitude 
read by a head varies Hneariy as the head moves off track 
center. Referring to Fig. 6(b), the read-back signal am- 
plitude from the head will be at its maximum during burst 
A when the head is positioned over track N + 2, N + 
4, • • • , As the head moves off track, the amplitude of 
the burst A signal decreases in proportion to the amount 
of otf-track position, and burst B amplitude begins to in- 
crease from 0. At tracks A^ + I . A^ + 3. A^ -I- 5, • - ■ , the 
signal from burst B is at its maximum and burst A is zero. 
Like the phase encoded servo, this scheme will tolerate 
small disk defects since many cycles arc used in the burst. 
However, amplitude servo is not as linear as a phase en- 
coded servo and there are two widely separated frequency 
components in the pattern— do and the burst frequency. 
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The H" nnodel matching algorithm Is proposed for the robust motor tracking control of 
linear time-varying multi-input, single-output uncertain system. This control method is 
verified by the system of the permanent magnet DC motor to improve the drive 
performance and reduce sensitivity of parameter variations, nonlinear effects, and other 
disturbances such as load changes and backlash. The matching model is used to identify 
the H"' tuning controller with bounded noise energy and minimized tracking error. An 
algorithm based on discrete Ricatti equation and H*" interpolation theory is used to solve 
the uncertain internal modeling and noise measurement issues in the tracking control 
system. The adaptive tuning controller effectively controls the presence of external 
disturbances, and identifies the rapid jumping and slow changing trajectory of the 
system. The effectiveness of the H« model matching controller for the motor tracking 
control system is demonstrated by using the developed adaptive tuning software 
controller. The H « model matching control method is examined through computer 
simulation, and then tested in a real time microprocessor-controlled drive system. 
Experimental results based on PC based microcomputer implementation are presented to 
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The Hinf Model Matching Approach to 
Motor Tracking Control System 



JeeMing Lee 
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ABSTR^ACT 

The model inaicliing algoriilim is proposed for 
the robust motor tracking control of linear linie-varying 
nuilti-input, single-output uncertain sysitin. This control 
method is verified by the system of the permanent 
niagnei DC motor to improve the diivc perforniajice and 
reduce sensitivity of parameter variations, nonlinear 
cfTects, and other disturbances sucii as load changes and 
backlash. I'he matching model is used to identify the H"^ 
Miming controller with bounded noise energy and 
mi^imi^ed tracking error An ali;orilhm based on 
discrete Kicaiii equation and 11^ inierpolauon ihcoiy is 
used to solve (lie nnceitain iiuernal inodelinii and noise 
meastirement issues in the tracking control system. Tlie 
adaptive mning controller effeci'ively controls the 
presence of external distnrbanc.es, and identifies the rapid 
jumping and slow changing trajectory of the sysieni. 
Tlie effectiveness of the model maicliing controller 
for the motor tracking control system is demonstrated by 
using the developed H'^-' adaptive (iming software 
controller. The model matching contiol method is 
examined through computer . simulation, and then rested 
in a real time microprocessor-controlled drive system. 
Experimental results based oit PC based niicrocompuier 
implementation are presented lo illustrate improved 
response and robust performance. 

Keywords: H"* model matching, tracking control 'system, 
bounded noise, adaptive tuning. 



I. INTRODUCT/ON 

Following the output of a system track from a given 
/efcreace trajectory is a common industrial problem. 



Sometimes the mechanical load vanes considerably, but 
the trajectory of the motor drive still has to follow the 
required tracking operation. These situations can be 
found in typical examples, such as: robot joints, and 
machine iools[ll[2]fS][9]f 10]. All these motion systems 
leqiiire the motor bore the load changing and still keep 
tlie prescribed trajectory. 

The motion disturbance in the above description 
accords to the nnccrtainiy noise under modeling system 
which has to bo minimized. It is the aim to make the 
motor drive system to meet the robust stability of control 
perfoimance. And in otder to obtain the satisfactory 
tiackingpcrformancel?), (he acbptivc tuning controller is 
usnally coni-idered in the linear on-line control stiaiegy 
for a :incccssful applicaiion( 1 1. However, it needs lo 
riusign a reasonably accuracy model, or the resulting 
control algorithn) is hard lo maintain the tracking system 
in the robust stability. Envisioned the performance 
exploited by robust stiibie modeling operation in the 
sense, it is believed that the controller should be 
designed to have abilities to bound the noise disturbance 
and mi-niniize the tracking error. 

In this piJ]>er, a motor tracking control design using 
H"^ model matching is proposed it is designed lo 
exclude the noise spectral densities whereas the 
conventional approaches can not. The merits of the 
approach he m the simplicity of the scheme and its 
practicality for real-time implementation. 

The identification has to be carried out within internal 
model matching condition. In order to make sure ihat a 
control transieni can be concluded already with the 
desired response, the identificaiion and ihe adaptation of 
the cofiirallers have to be peifonned as fast as possible. 
The limited time needed for updating the coiitrollers is 
determined by the convergence time of Ihe identification 
algoritiiin and the lime needed to calculate the new 
controller sutting. 
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Generally, clcctricn! drives are nonlinear Kj'stems 
because of current, vollage and speed boundaries. In 
addition, their parameters, e.g. resistance, inductance, 
inomeiii of inerlia. can vary. Also in real drive systems, 
the disturbances, e.g. load torqiie on the motor sliaft. 
must be considered The non-linearity of load torque 
and njechanical dampers is the main factors to make the 
tracking misalignment and system unsrable[7J. In this 
paper, the system model of pcrnianent-magnet DC motor 
is adopted to verify the tracking performance and robust 
stabih'ty. 

Fig.l-I is the block scheme diagram of main sections. 
Baishlcss dc motor represents in the dotted box. 
Inverted conunutation circiiit shows uitlnn the dashed 
line. The others of the diagram are included with the 
controller. Two input signals of the system from the 
cnrrenl sensittg which is proportional to the motor 
voUage and (he increment encoder determine the rotor 
posidon 
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stator becomes zero, making the dynamic model simpler 
The relationship between voltage and current of a 
permanent magnet motor drive is as follows: 



n 



Usually, iq make the minor cfTcct in the controlled motor 
system. Therefore, the id current of the motor produce 
the strong control factor in the system. From eqn (2-0 
Ed, Eq arc the motor's terminal voltage in the steady 
state. The dynamic electrical circuit model is derived 
as: 

d i^, fill - -RnlLa.^j - (5 (OjU + E„/A„ (2-2) 
The instanlnneoiis torque '// is 

7) - Ji (bm ^- + k(f^^ (2-3) 
The system parameter is described as follows: 

/?„ ; the armature resistance 

L„ : the armature inductance 

Ji : ihe motor moment of inerlia 

^J; the magnetic flux of the permanent magnet 

(constant), 
w/w : Ihe angular velocity of the motor 
Bfti : the coefficient of viscous friction 

Ti : the load torque 

Ed, Eq, i(f and iq are d,q shafl voltages and currents. 

From eqn (2-2) and (2-3) (he transfer function of Ihe 
permanent magnet motor drive with the variable 
disturbance load is shown in Fig.2-1. 

The dynamic model of pennanent magnet DC motor is 
well established in the continuous time domain by Ihe 
state space modet[13]: 

x(l) = Ax(l) + BU(t) (2-4) 
where 



BRUSHLESS DC MOTOR-: 



x(t)=[e„,(0 a)„(/) J AO igin] 

U(t) = [n(0 BAO E^(0] 



Pig. 1-1 The hnrdwme scheme diagram of motor main 
system 
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n. DYNAMIC MODEL OF PERMANENT MAGNET 
MOTOR DRIVE SYSTEM 

The permanent magnet motor drive can be regarded 
as direct current by having d-q axes transformation that 
rotates synchronously with the alternating current. 
Because the relative velocity l5ctv\'een the rotor sliaft and 
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where 

x(t) = state vector of the motor components. 
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U(t) = the disturbance (load torque) and the control 
signal (voitags) 

A and B are sysleni niaf rices corresponding to the 
iiipuL/oiifpiit veclors. respectively. 

''1 I 



it 



Fig. 2- 1 Block diagram of a bnjsliless DC motor 

To protect tlie niolor from overload in t)ic motor 
tracking system. Tlie niolor system has specified 
constrainis on the motor current. These limitations are 
a.<;foIlo\vsf2Jf71[Ill: 

(1) Peak current limilation: The niaxiniiini occurring 
ntolor is liniiled (o avoid demagneti7.ation of the 
permanent magnets. 

(2) CoiUiniious current limitation To prevent the nmtor 
from overheating, tlie maxinuini constant current is 
limited to appro.ximateW half the peak current for the 
motor used. The continuous current limitation becomes 
active only after a certain period of time and can be 
disregarded in position control. However, if torque 
control is used, longer jxiriods of maximum applied 
power can occur; in this case, tlie continuous current 
limited is important for the prevention of the power 
capacity of (he driver coniponents. 

f!J) ("onmiutation limitation: The inrush current through 
the power drive components increases with the motor 
.ipeed and (he heavy load varialion. To prevent the molor 
from being damaged, the product of armature current and 
voltage is limited. 

In this paper, however, a systematic method is 
proposed to determine the plant orders and time delay of 
(he motor tracking systemf I ][1 7]. To derive the 
computation equations, (he time domain model of PM 
machine with a controller is transformed into the (he 
discrete model in the z domain form, and generates the 
system upper bound for tlie worse case, as well as robust 
tracking output. The system is specified the requirement 
as foIloHings: 

speed bandwidth: 0 - .5OH7, , imder load condition and 
l2.5Hz for loading variation. The robustness is as large 
as possible within (lie setting frequency response. The 
accuracy of the position is within lum, and no exciting 
phenomenon. Inscnsitivily of load disturbance and noisy 



interference is designed in (he adaptive control scheme 
and modeling inatching. 



in. THE H*^ SYNTHESIS CONTROL CONCEPTS 

Firs(. the H*^ syn(hesis control system is structured as a 
standard problem, depicted in Figure 3-1. where the 
objective is to design a real-rational proper controller K 
such that tlie closed-loop system formed by real rational 
proper G and K is internally stable and the norm of 
the transfer function between external input vector v and 
controlled output vector ? is minimized over all 
stabilizing controllersf3]. 



Fig. 3-1 H"" synthesis standard problem 



H' : exogenous input vector included disturbance d, 
noise n. and reference order r. 

2 : controlled output vector included tracking error 
Te. energy z^, and output tracking z^ 

y : observed measured output vector. 

u ; controlled input State. 

G : The augmented control model. 

K : The feedback controller 

This augmented control system C can be written as: 



G: r(<7/-/l) *^-f D 

The synthesis system can be represented as: 

Setting tlie system gain K for feedback as: 

V " Ky 



(3-n 



f3-2) 



(3-3) 

The convergence and stability of the system can be 
obtained from inputs w to the outputs z of tlie close loop 
function. The form is expressed as: 

z-'-TnAO.nw (3-4) 

^ (Ju + Gi2K(/-a2iK)-^G2i (3-5) 
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Eqn(3-5) represeiils the standard form of H'^ synthesis 
control system. Tlie augmented synthesis s>'sletn has to 
find the controDer K so that 

inf||cMi \ nnk'iJ -GyifCyUlu || (3-6) 

The controller K [4J has to satisfy the following tvvo 
conditions: 

(i) The closed loop system is internally stable. 

(ii) minj|7-„r((';.^')|L <y 0<v<| 

The augmented control system G means the combination 
of (he real control s>'Stem P and Hltcr functions aW for 
the tracking error and the internal system entropy. The 
filler fimction wc adopt in this article shall consider llie 
response of the oiitptit tracking as well as the minimized 
sensitivity for the system stability thereby boosting the 
system's performance and robustness. 

The synthe.sis Tl"^ system is examined in 
reconstruction slate space forms[5)f6]. The mean square 
reconstr\tction error value is to measure Die observer's 
reconslniction status. To design a time varying tuning 
procedure as follows: 

Qin I I)- /I(/0(>'(/0y(/n + /f',{'7) 

( p re- CO va r i a n cc ) (3-7) 

l.{n-^\)- <J{n + {n ^ \)[CXn -t- l)(;(n)r' + 1) 
t ' (state estimation gain) (3-8) 

(observed covari an ce) (3-9) 

From Ihe above adaptive procedure, the synthesis H'^ 
system shall be examined from (he predictor-corrector 
state. 

x{n t I)" A(n)xiit)^- fHn)u{n) ^ Lin A !){ 

yiu V 1) •nM)f/l(j7).?(;7)+i?{7r)H(/0]} (3-IO^» 



IV. THE ADAPTIVE MODELING IDENT mCATION 

Tiic dynamic model of PM DC motor is ^^'ell set up by 
the augmented system that is translated into the 
difference equation necessary for use in the discrete 
systemf 15]. 

The discreic form is given by eqn (4-1) in the state 
space model: 

V, ^ A in - \ )x(n ■{)+ 7f{fi - [W(n - I) (4- !) 

The matrices /I fn- I), and /?(;7- 1) in discrete-time 
model are based on the motor's parameters in cqn (2-4). 

The approach to make the identification of a motor 
trajectory employs the QR based method to solve tlie 
recursive least square(RLS) problem. , The QR based 
tnethod is numerically stable. It shows that the 



computation mclhod of matrix problems is for orthogonal 
similarity transformation. 

To do so in the adaptive tuning proccdurel 1 7], the 
objective is to fmally obtain a predictor model for the 
desired tracking v, :. 

y,(n) = X'^Oi)B{n) + C, (4-1) 

where X'^(t) is an observed vector using past values of 
yi and X. BCt) is a parameter vector using the coefficients 
from the above of A(n-{) and A?(h-I), and Q is an error 
variable representing the load disturbance and modeling 
error. From the based least squares computation, it 
is to see that generation of pseudo-linear variable Vi. 
defined as 

y^^r'ATXCn) (4-2) 

The setting up recursive least square (RLS) model, keep 
the identification order low. A state space predictor for 
the plant modelfMJ can now be set up by a series of 
estimation steps. 

error t^, = - .V(»)9('0 (4-3) 

controller gain; 

/^;(") - 'In /;.,oo^+ \„,(")r' /;_,(n);i 

(4^4) 

estimate: e^(n) = e,_,(»)+ Al'^_,(/Oc, (4-5) 
forgetting factor: 

X, («) -/-(/ + fY (« / r (4-6) 

covariance: 

/;(n) = M I />A://or/;_,(/r)f.U^A',(w)i 

^^/(")^(")/^,(") (4-7) 

/',(«) has to be positive semidefinile. and \Vi*,{n) has rank 
defect r. then it needs to use an orthogonal 
transformation Tto adjust it. 

V. The Model Matching for llie Motor System 

1 he estimator of Ihe system is reconstructed after 
every periodic computation. It consists of a set of time 
varying equations from eEqn (3-7) to eqn (3-10). These 
works can form an estimated system model. The 
controller of the augmented system is designed in 
accordance with eqn (4-2) to eqn (4-7). It executes in 
Ihe state of the adaptive tuning to minimize the motor 
tracking error. 

The approach design is needed to specify the 
disturbance from the input for the worse case. Therefore, 
the term ganmia y is required to setup for the robust 
range. 

<y{hOOQuo)<y^ (5-0 
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Fig. 5-1: The synthesis model control .Tysteni 

It means Ilial the spectral radius of P,{n)Q(n) must be 
lower than the upper bound of disturbance V- 

Fig.5-1 shows the synthesis control scheme model 
for the motor control system. From the figure 5-1 the 
estimated system model is derived from the real motor 
system and the desired augmented reference plant. Tlic 
adaptive controller is desigiied after the system 
model 116] is estimated. 



dotted line on Fig. 6-1 and Fig. 6-2 represdtlts the deilrtil 
trajectory, while the the dash line represeiili iht mtpiii 
trajectory. The position trajectory follows Ihft d^lred 
signioidal tracking on Fig. 6-1. The speed trajectory 
gets some disturbance but still maintain the set speed 



200 

1 



Fig. 6-1 Case I: Position trajectory 



VI. TFuST RESULTS 



Fig. 6-2 Case 1: Speed trajectory 



Several test cases in (he experiment evaluate ihc 
proposed H"" approach model matching technique. 

The position track selected in lhis test is composed of 
three segments of si^moidal functions. The first rotates 
the fonv^rd position by 200 radian, the second is holding 
the position at 200 radian for 2 second, and the third 
turns the reverse direction back to the original position 
by the armature current braking from (he feedback sensor 
signal. In these test cases, the armature surge current is 
limited within the specified range to verify the motor 
system nmning uithin the designed trackijig and safety 
working condition 

Two cases are shown on the following mainly to 
represent the trajectories of the motor rotor position and 
requirrid speed tracking paths. Each of the motor system 
tracking is measured from (he starting time -and sampled 
in every 20 msec. 

Case f shows figures from Fig. 6- 1 to Fig. 6- 10 to examine 
(hp motor system performance and robust. Case 1 
represents 10% perturbation with the required trajectory 
entering the ntotor system and the load torque changing 
periodically every I2.5Hz, Figure 6-1 to 6-7 shows the 
tracks and the measurements of the rotor position, speed, 
voltage, and electric torque, position tolerance, direct 
current, and quardrature current, respectively. The 
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Fig. 6-3 Case 1: Position tolerance 




Fig. 6-4 Case I: Electric torque 
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Figure 6-5 Case I: Armature voltage 
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Fig. 6-9 Case I: Parameters of B(q) 
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Fig. 6-6 Case I: Direct current 
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Fig, 6-7 Case I: Quardratu re current 
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Fig. 6-8 Case I: Parameters of A(q) 
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Fig r>-10 Case I: Parameter of C(q) 

running on Fig.6-2. Fig. 6-4 shows the electric torque 
increased during the motor accelerating (he speed. 
However, the motor system stiti maintains in the safe 
running condition by examining the voltage and current 
variation from Fig. 6-5 to Fig. 6-7. 

The time varying motor system under the noise 
perturbation and load change highly depends on the 
adaptive tuning for the satble and robust synthesis 
system performance. Fig. 6-8 to Fig. 6-10 show the 
synthesis parameters of the system. The parameters 
of matrices of A(qX B(q) and C(q) show the adaptive 
identification for the reduction order H"* approach 
matching. Specifically, when the motor proceeds the 
tracking in the braking segment, the system parameters 
vary quite largely. 

To satify the robust performance of the H" synthesis 
approach, it is designed the upper bound for the worse 
case at the gamma value (y) 0.67 in Case F. 

Case IJ represents the situation of the 50% uniform 
distributed perturbation with the desired trajectory 
entering the motor system and the load torque changing 
periodically every 25Hz. Fig. 6-11 to Fig. 6-15 shows 
the performance of the motor system in case 11. It 
examines tlie motor system with tlie robust control in the 
H«' synthesis approach. In Case IK the gamma term y is 
designed 0.45 for the robust range. It can meet the case 
to limit the uppier bound of disturbance. 
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Figure 6- M Case IP Position trajectory 



Fig. 6-12 Case II: Speed trajectory 
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Fig. 6-15 Case II: Electric torque 

noise and load \'arying. It can be seen from Pig. 6-13 
tliat the position trajectory has one step predicted error, 
especially, during the motor increasing speed. However, 
the dc error of the position trajectory is derived from 
the uncertainty disturbance. It can be compensated to 
adjust the deadbeat of the filter 

In this paper, it mainly empliasizes the performance 
of the position trajectory. Therefore, the speed 
Irajcctoryon Fig. 6-12 is infected by some disturtence 
siniilarly as Case I. during the motor starling, speedup, 
and braking stop. Fig. 14 and Fig. 15 show the motor 
system still stable when the strong perturbation adds to 
the system In Case II it stiil can work within the 
protected safty range. 

Fig. 0-13 shows that the larger position tolerance appears 
as the disturbance become larger, comparing with Case 
1. It shows that the dc tolerance is reduced by adding the 
constant backlash. However, it is designed from I he 
."Ty'nthesis s)'stem to bound the disturbance in the spectral 
radius. 



Figure 6- 1 3 Case li: Position tolerance 



CONCLUSIONS 



Fig. 0-M Case II: Armature voltage 

The dolled line on Fig 6-11 and Fig. 6-12 represents 
the desired trajectory, while the the dbsh line represents 
the oulput trajectory. Fig 6-11 shows the position 
trajectory following the desired tracking under the heavy 



In this paper, the design and implementation of an H"' 
matching controller with an adaptive tuning has been 
investigated A syntiiesis systematic and analysis 
approach for developing the algorithm is presented. The 
H"^ filter is tuned by various parameters 
Qi")JAr),K(f7)M").Pi(nlofid X(n) which are 
dominant in the slate identified model and the designed 
H'*' controller. These paramelers have to meet the 
dciired transient and stead>' slate pcrfbrrnance. The 
discrete adaptive H"^ matching controller has been 
designed to be well suited the motor speed and position 
trajectory for Bnishless DC motor. The load torque 
variation and mechanical damper can be immunize by 
the designed H'^ controller under the upper bound area. 
The ne.xt step in this research will lest the FT* matching 



- 630 - 



conl roller in difTe.re.nt kinds of disturbance Tor tlic speed 
nnd posilion trajectory. 
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Nonlinear State Observation Using 
ifcxj-Filtering Riccati Design 

Konrad Rcif, Frank Sonnemana and Rolf Unbehauen 



Abstract — In this paper, the authors propose an observer for 
continuous-rime nonJincar .syjfteras. The observer gain is computed by 
a Riccati differential equation similar to the extended Kalman filter. 
They prove that under certain conditions the proposed obsener is mi 
cxponeiirial observer by choosing an appropriate Lyapunov function. 
Furthermore, the authors explore some important reiations of the 
proposed observer to robust control thcary and Hoc -filtering. To 
examine the practical usefulness of the proposed observer they applied 
It to an induction motor for the estimation of the rotor flux and the 
angular velocity. 

Index Terms — Extended Kalman filter* ifo^ -filtering, nonlinear oIh 
ser\cr. 



I. Introduction 

The basic model for the general nonlinear observation problem 
consists of a nonlinear plant and a nonlinear output map, at which 
neither system noise nor observation noise is considered. The problem 
is to estimate the states of the plant by observing the output when 
the initial state of the plant is not known. In literature the design 
of nonhnear observers is a widespread area of current research (see 
[13] and [22] for a survey). Altliough many important results have 
been achieved, there are still possibilities for an improvement of the 
existing synthesis methods, 

A frequently apphed method is the extended Kalman filter [5], used 
as a deterministic observer for nonlinear systems. The advantage is 
its easy handling with a simultaneous good performance; moreover, 
a full proof of its exponential stability has been given recently [16], 
[1 8]. However, a well-known difficulty arising through its application 
is the problem of divergence. The reason is that in computation of 
the observer gain only the linearized funclioni; are considered^ but 
not the whole nonlinearities. To obtain convergence, it is therefore, 
in general, necessary either to have a good hiitial guess so that the 
initial estimation error is sufQciently small or to permit functions 
only weaklN*^ nonlinear [16], [ItS]. 

Recently ffoc -control methods have attracted attention in robust 
control theory [4], [1 ij. It has been mentioned that the high robustness 
of tlie LQ-optimal control systems can be further improved by using 
^cc: -control [3], In this paper we use a .similar approach to improve 
the extended Kalman filter. Motivated by the results of robust control 
theory [9], [20] and /^cv^-filtering [2], [6], [10], [14]. The remaining 
nonlinear terms are treated as uncertainties and are implemented in 
tlie proposed observer design by a worst case valuation in such a 
way that tliese uncertainties can be tolerated. Consequently, better 
convergence properties can be expected; moreover, under certain 
conditions convergence can be achieved even independent of the 
initial estimation error. An alternative method for improving the 

.Manuscript received March 25, 1990; revised September 23, 1996 and 
March 18, 1997. This work was supported in part by the Univcrsitut Erlaiigcn- 
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convergence properties of tlie extended Kalman filter has been 
proposed in [17], 

The paper is organized as follows: in Section II we state the nec- 
essary preliminaries and introduce the proposed nonlinear observer. 
By choosing an appropriate U^ipunov function we .show that the 
proposed observer is an exponential observer in the sense of [8]. 
Important relations to if cc -filtering [2], [6], [10], [14] are explored in 
Section III. In Section {V we apply the proposed obser\^er to estimate 
the rotor flux and the angular velocity of an induction motor (cf, [12], 
[15], [19]). Some conclusions are drawn in Section V. 

11. Proposed Observer 
Consider a nonlinear system represented by the equations 

m = cz{t) (2) 

where z{t) € FV^ is die state. x{t) G the input, y{t) G i?."' the 
output, and C a m x q matrix. The function and the input 

x(-) are assmncd to be C^-fmictions. For the sake of simplicity we 
arc restricted to linear time-invariant readout maps. The important 
case where some of the .state wiablcs can be measured directly is 
thereby treated. Furtiiermore, tltc linear time-invariant readout map is 
well suited to pointing out the main idea. Some remarks on the more 
general case with nonlinear readout maps arc given in Section HI. For 
the considerations following below we need a slightly more general 
definition of the exponential observer in [8]. 

Definition: For a dynamic system given by (1) and (2) an expo- 
nential observer is a dynamic system 

i(0 = f(m-Mt)) + 9W)r'(:)M-ht) (3) 
with an output injection <;(•,•. •,•) such that 

\m-m\\<-n\\m-m¥-'' 

Vf > 0, V^(0) - f(a) £ (4) 

holds for some p,i],6 G with Bf, = {v G R'^ : \\v\\ < p]. 

Remark: The output injection g{'. •) depends in general on the 
whole functions i(') and x(') and not only on the recent values z{t) 
and Tlie reasons are given below. 

To calculate the output injection •) we introduce a Riccati 

differential equation which is similar to the Riccati differential 
equation for the extended Kalman filter [5]. Using the abbreviation 

A(0 = |^(i(/)..rU)) (5) 

it reads 

P{t} = A{t)P{i)^P{i).e{t) 

- P{t){C^R-'C - A-i)P(f) -h fi-r (6) 

with some A, // 6 R^ and a positive definite symmetric matrix R. 
We define tlie observer gain by 

K{t.)^P{t)C'^^R'' (7) 

and the output injection by^ 

i/{yit),M').xi-)A) = /mOMO- C£(f,)). (8) 

It is necessary to allow a general dependence on and x{') and 
not only on the recent values z(i) and x{t) because the solution 
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of the Riccati differential equation (6) depends on past values of 
.4(0 = (df/ds){m^x(t)). 
An expansion of f(\ ■) at z{ t), x{t) leads to 

= A{t){z(t} - i(0) + ^(::(t),z(t),x(t)) (9) 

where ipU{tXz{t),:i-{t)) are the temis of second- and higher order 
in - i(0. Defining the esiimation error by 

CfO = ^(0 " m (10) 

subtracting (3) from (1) and considering (8) and (9), we obtain 

at) = {A{t) - K{t)C)at)-\~9(^itlz(t},:r{t)). (11) 

For an exponential obsen^er tliis differential equation has an expo- 
neniially stable equilibrium point at zero; moreover, if (4) holds for 
all C(0) 6 i?^ the equilibrium point is globally exponentially stable 
(cf. [21, Sec. 5.1, p. 143]). 

To state the proof of the Main theorem given below we need the 
following preparation. 

Lamina: Consider real vectors z,z G R"^ and ;r € a symmetric 
q X q matrix 11, and a function with 



(12) 



for some k ^ R^. Choose A,/i- E R ' such that k. < A;/ holds. Then 
there is a i > 0 such that 

2{z - zfn-Az, Lx) <{z- zf[x'i + (/r - ^)nn](.- - i) 

(13) 

holds. 

Proof: From k < A// it follows that there is a 0 < ^' < A// 
such (hat 



t\ < A// — <!>" 
holds. Inserting inio (12) yields 

Mz,z.T)\\<{X,,-6')\\z-i\\ 
and therefore we have 

■2{- - ifn^(z,i..r) < 2{\m - *')l|n(3 - - 
For considerations following below we evaluate the square 



(14) 



(J^,-Ly\nu-B)\\-x\\:-2\\j 

r r. 

2//* . 6 



- 2A (^/' - ^ j \mz - z)\\\U - i\\ + A^^- - i)''!-- - i) > 0 

and use the fact that it takes only nonnegative values. Rearranging 
the terms yields 



2(V-5)||n(z-£;||i|(^-i) 



Choosing 



nn + A'i 



(15) 



6 = 



2n6 



A^ 



from 8 < \fi- it follows that fj > S /\ md therefore 
6 fi 6- 

Combining (14) and (15) we obtain the de^sired result (13). □ 
With these preliminaries we are able to state Theorem 1 . 
Theo7'em J: Consider a nonlinear system witii linear rime-invariant 
readout map given by (1 ), (2) and an observer given by (3) and (6)-(8) 
with some A. ft > 0 and a positive definite matrix 7?. such that the 
following condirions hold. 

1) For any solution i(-) of the obsen^er differential equation (3) 
the solurion P(-) of die Riccati differential equation (6) is 
bounded via 



pi < P(t) < pi 



(16) 



for some p,p > 0. 
2} The nonlinear terms ^{z{t), z(t).x{t)) in the difTcrential equa- 
tion (1 1) for the estimarion error are bounded via 



||v.(^(f), £(0.40)11 <^ll--(^)- i(OII 



(17) 



with < A/f. 

Then the considered observer is an exponential observ'er. Moreover, 
the differential equation (11) for tlie estimation error is globally 
exponentially stable at zero. 

Remarks: 

1 ) To ensure the global exponential stability of the differential 
equarion (11) we require tliat (17) is satisfied in the whole 
stale space. If this inequality holds oniy in a restricted area, only 
local exponential stability can be guaranteed. This means that in 
this case tlie initial estimation error may not be arbitrary large. 
Nevertheless, numerical simulations have shown that even if 
the stability is not globally satisfied, the proposed observer can 
tolerate much higher initial estimation errors than an observer 
based on the extended Kalmaji filter (see Section IV). 

2) Using similar techniques as in the proof following below, it 
can be shown that under certain conditions the continuous-time 
extended Kalman filter is an exponential observer [16], [18] too. 
The difference is that the equilibrium point of the estimation 
eiTor is globally exponentially stable for the proposed observer 
but not for die extended Kalman filter. Tliis means that the 
proposed observer can tolerate arbitrary large initial estimation 
errors, which is not the case for the extended Kahnan filter. 

3) For tlie usual extended Kaimaii filter it can be shown that the 
solution of the Riccati differential equation is bounded via (16), 
if the matrices A{i) and C fulfill the uniform observability 
condition [7] or if the nonlinear system satisfies a certain 
deiectability condition as stated in [1]. Unfortunately, these 
conditions cannot be applied to get similar results for the 
ifo:; -filtering-like Riccati differential equation (6). (The reason 
is that tiie term R~^C — A"/ is in general not positive 
definite because R~^C docs not have full rank if we n-eai 
nontrivial observ^ation problems.) It is a well-known problem 
arising in H ^-control as well as if^c -filtering that tlie solutions 
of the corresponding algebraic Riccati equations and Riccati 
differential equarions may lack being posirive definite (cf. [3]). 

4) For the proof below as well as for the lemma stated above, 
the parameter does not play any important role. It has been 
introduced to have an addirionaJ degree of freedom during the 
implementation process. Especially for ensuring Condition J, 
the additional design parameter //, has turned out to be very 
useful. Also see tlie following remark, 

5) The existence of a solution P(-) for the Riccati differential 
equation (6) which satisfies Condition 1 depends mainly on an 
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Using 



K{t)C(f.) 



A{t) 



<fiz(i).S{t).x(t)) 



Fig. 1. Block diagram corresponding to the diffcrcnliiil equation (27). 
v,it) 



[m] 
[ml 



v(t.) 



Id -K{1)]^ 



A{t)-K{i)Cit) 



linear part 



Fig. 2. Block diagram corresponding to the difTereiitial equation (31). The 
whole system consists of two interconnected sub.wstems. Tfic dashed box 
contains the linear part of tfjc system and the observer; Fit) and Hit) 
represent the higher order nonlinear terms. 



appropriate choice of A and jl As will be shown in Section III, 
they arc closely related lo the attenuation con.stant 7 in H^- 
filiering problems. To find suitable values for A and }t we 
employed a binary search algorithm, which is widely used to 
solve i^oc -control and Hcx.-iiitering prcblenis (see, e.g., [3]) 
and adapted it for our purpose (see Section IV). 
Proof: We consider the differential equation (11) for the es- 
timation error and prove its exponential stability by choosing the 
Lyapunov function 



v'(C(0.0-c'V)n(z)c(.0 



(18) 



with nu) = P ^(t). Because of Assumption 1 we have the 
following bounds for tlie Lyapunov function: 



M)f <ymA)<-\\mf 

V V 



(19) 



which imply positive defmiteness and decrcscence. The time deriva- 
tive of the Lyapunov fiuiction reads 

Inserting (,'(/) according to the differential equation (11) yields 
Applying the lemma and considering (7) we get with k < A/i 

- 2c''r-'c + (//' - 6)n(on{o + x'm.t)- 



n(0 = ~E{t)P{t)U(t) 



(20) 



and the Riccati differential equation (6) leads to 

y{C{tU) < C''\t.)[-6n(t)n{t.) - C'^'R-'Qat), (21) 
With the bounds (16) for P(t) we obtain the inequality 



(22) 



which implies that the time derivative of the Lyapunov function is 
negative definite. With (19) and (22) we satisfy the requirements to 
apply [21, Th. 62, Sec. 5.3, p. 173] and conclude that the differential 
equation (11) has a globally exponentially stable equilibrium point 
at zero. This includes tliat the considered observer is an exponential 
observer. □ 

III. Relations to -Filtering and Robust Control Theory 
Considering the small gain liieorein the proposed obser\'er is 
closely related to IToo-filtering theory as will be sho\^Ti below. Since 
these relations are more obvious for the case with nonlinear readout 
maps, we let a nonlinear continuous-time system be given by 

i{t) = f{z(t),x(f)) (23) 
V(t) = (24) 

with a -function /j.('). The linear readout map is thereby contained 
as a special case. For this nonlinear system we construct an observer 
according to (3) and (6)-(8). Using the fact that tlie functions /(•, •) 
and //(-) are continuously differentiable, we expand them via (9) and 

k{z{t)) - hizit)) = c(t)(z{t) - ~{t)) -h x(m.m) (25) 

where C{i) is an m x q matrLx given by 



c-(0 = |^(i(0). 



(26) 



Defining the estimalion error by (10) we obtain from (23) and (3) 
using (8). (9), and (25) 

C(*) = (A(t) - K[t)C(t))C{t) + 9(3(0. i(0,^-(0) 

- K(t)x[zit)^Mt)). (27) 

The nonlinear terms are assumed to satisfy (17) as well as 

\\x(m-Htm<\\Ht)-z(t)\\ (28) 

with some A > 0. Now we have a closer look at the differential 
equation (27) for the estimation error. Fig. 1 shows the corresponding 
block diagram. Defining 

^ -j9(5(0,i(t),-T(f))c''(f; 



F{f) = 



/'IK(.<)II^ 



(29) 
(30) 



mm 

for at) ^ 0 we have from (27) 

00 = {A(t) - K{t)C(t) + DF(t) - K{t)H{f)K(t) (31) 
with D = fil as shown in Fig. 2. Moreover, wc obtain 

ll[S'.!]||<^ 

for every t > Q which is an immediate consequence of (17), 
(28)-<30), and h < A/t. By this construction the nonlinearities are 
considered as time-vary^ing uncertainties^. For the stability of (31), we 
have the following theorem. 
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Theorem 2: Consider (31) with D = /// and an observer gain 
Kit) given by (6) and (7) with A,/i > 0, and R - 1. Let the 
following aindilions hold. 

1) The Riccati differential equation (6) has a solution P(-) which 
is bounded via 

pi < ^'ij) < pi 

for some j?./7 > U, 

2) The matrices F{t) and H{t) are norm-bounded by 



Fit) 
Hit) J 



< A 



for every t > 0. 
Then tlie system (31) is exponentially stable. 

Proof: For the proof of ihis theorem we employ a variant of 
the small gain theorem. We recall the operator nomis [10] for the 
uvo subsystems in Fig. 2 (e.g., [F^ (t) ff '^ [t)]^ and tlie pai1 inside 
the dashed box) 



||ru.v||.:x3 = sup 



||r...||.^= sup PI 



where || - II2 denotes the usual Lo-norm for functions. According to 
[20, Corollary 1], the whole system in Fig. 2 is exponentially stable 
if the subsystem in the dashed box is exponentially stable and the 
small gain condition 

||r...||oc(|7;,,||.^ < 1 (32) 
is satisfied. The operator nonn ||r^.v||oe can be estimated easily to 

(33) 



irru>v||.aj = sup 
f>0 



m\\\ 



Moreover, it c^n be seen that the subsystem inside the da,shed box 
in Fig. 2 has the structure of a -filter. Introducing 7 = 1/A and 



using R — I, D — It I the Riccati differential equation (6) can be 
rewritten in tlie form 

Pii) = A{l)P{t)^~P{t)A'it) 

- Pit)C"^'(t)C{t)P(t) + + DD'^ (34) 



and die observer gain (7) is given by 

Kit) = P(f)C^'((), 



(35) 



With (34) and (35) we are able to employ standard results in H^- 
filtering (cf. [2, Sec. 7.4, p. 301, Th. 7.6] or [6, Sec. 7.2.2, p. 267, 
Th. 7.2.1]) and conclude that the subsystem inside the dashed box 
in Fig. 2 is exponentially stable and tlie operator norm ||rru,'|iw is 
bounded by 

lir,.u.|U. < 7. 

Together with (33) and 7 = 1/A vv'e fulfill the small gain condition 
(32) and the proof is completed □ 

IV. APPLICATION TO AN EXAMPLE SYSTEM 

To examine the practical usefulness of tlic proposed observer we 
apply it to a liighiy nonlinear example system, a symmetrical three- 
phase induction motor. The quantities to be estimated are the rotor 
flux and the angular velocity (see also [19] znd the references cited 
therem). We introduce the state differential equations for an induction 
machine and do not mention too many details about the physical 
meanings of the considered quantities, which are contained in any 
standard book about electrical diives (see, e.g., [12, Sec. 10,1] or 
[15, Sec. 3.3]). The .slate differential equations of a symmetrical 
three-phase induction machine read 

ilit) = kiZi{t)~^Xi{t)Z2{t}-^k2Z3{t]-\-X2{t) (36) 

hit) = ~xj{t)z,it) + hzuit) -f h2Z.^it) (37) 
= k-3Ziit) -H h-iZiit) + (;ri(^) - z,,{t))z.iit) (38) 



IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 44, NO. L JANUARY 1999 



207 




10 15 20 

normalised time 



25 



30 



Fig. 4. Lyapunov function V(X(t),t) for the proposed obset-ver. 



hit) = k-sZ-At) ~ ixiit) - Mt)>At) + k;z4t) (39) 

i-,(t) = ko(zi{t)zAt)- Z2{t)z-4t))+ kiSX3[t). (40) 

All state variables are nonnalized ^i(0' ^^(t') and z-^it), z.i(t} are 
ihe components of ihe staior and the rotor flux, respecliveiy, in the 
plane perpendicular to the rotation axis. Zfj{i) is the angular velocity, 
xi(f) the frequency, and ;i'2(f) the amplitude of the stator voltage, 
respectively, and xiiit) is the torque of the load. A;i , ■ • • , are several 
parameters depending on the considered drive. 

The output are the normalized stator currents yi{t) and 1/2(0' 
leading to the output equations 



V2(t)^k7.Z2{t)-tk^z.i{t) 



(41) 
(42) 



\sath parameters A-7 and ka. The matrix A{t.) can be calculated from 
(36>-(40), parameter t omitted 



- k, 




k2 


0 


0 1 


~:ci 


ki 


0 


k2 


0 


k3 


0 


Am 


xi - 35 


— -4 


0 


k-. 


h - J-i 








— k^ss 


-'ks.^2 


kszi 


0 - 



A{t) : 



and the output matrix is given by (41), (42) 



C = 



0 

k7 



0 



(43) 



(44) 



For the simulation we set ki = -0.186, A--> = 0.178, ki = 0,225, 
Ah = -0.234. kr, = -O.OSL ke = -O.OIS, kr = 4.643, k^ = 
-4.44S, xiit) ^ .v(0. X2{f) = s{t), and x:s{t) = 0, where is 
the unit step function. These values are obtained, if in the original 
unnormaiized system ceilain values for ilie physical quantities are 
chosen. To obtain 'tq:propriate values for A and fi we propose a binary 
search algoritlim consisting of the following steps. 
1) Choose A,// > 0 and = rl with some r > 0. 



2) Decrease by taking the half-value until the solution P(-) 
of the Riccati differential equation (6) takes positive definite 
values. If no positive definite solution can be achieved, proceed 
Hdth Step 3). 

3) Decrease /t by taking the half- value and initialize R = rl 
as in Step 1). Proceed with Step 2). If gets smaller than a 
prescribed bound /^o we let the algorithm abort. 

For the simulaiioru; we insert A{t) and C into the Riccati differen- 
tial equation (6) and choose =r 0,7, /t^ = O.OOS, and R = 0.06 J 
aceording to tlie above stated search algorithm. The observer gain 
is given by (7) and tlie differential equation for the observer by (3) 
and (8). We .solved the state differential equations numerically by the 
Runge-Kutta metliod choosing the initial conditions 



.^(0) = [0.2 -0.6 -0.4 0-1 0-3]'^ 
for tlie system to be obser\'ed and 

z{{)) = [0.5 0.1 0.3 -0.2 4f 



(45) 



(46) 



for the observ^er. For the sake of comparison we also computed a usual 
extended Kalman filter [5] with the .same weighting matrices /rl, R 
and the same initial condition i(0) by .setting A = 0. fr = 0.008, 
and R = 0.06/ in (6). 

Figs. 3 and 4 show die simulation results. As shown in Fig. 3 the 
estimaiion error approaches zero for the propo.sed observer but not 
for extended Kalman filter. Further simulation results have shown 
that ji.?' proposed observer has often a larger domain of attraction than 
the extended Kalman filter, i.e., the initial estimation error may be 
larger. The Lyapunov-function (18) is depicted in Fig. 4. Looking 
at the values for 0 < /• < I we see that tiie Lyapunov function is 
increasing. This means that (22) is .sufficient but not neces.sary for 
the en-or decay. 
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V. Conclusion 
Tlie purpose of this paper is to propose an observer for nonlinear 
systems. This is done by modifying (lie well-known extended Kalman 
filter [5] for an effective treaaneut of the nonlinearities. We proved 
that under certain conditions this observer is an exponential observer 
[8] by choosing an appropriate Lyapimov function. For an alternative 
approach to the proposed observer we consider the nonlinearities 
as uncertainties and get a stabilization problem which is typical 
in robust control Uieory [9], [20]. This stabihzation problem can 
be solved by i^oc- -filtering techniques [2], [6J, [10], [14] leading 
to the proposed modification of tlie extended Kaiman filter. Using 
the proposed observer to estimate the rotor flux and the angular 
velocity of an induction motor we show that it can be easily applied 
to complex highly nonlinear systems. Numerical simulations show 
a good performance and an increased domain of convergence in 
comparison to the extended Kahuan filter, i.e., we can tolerate 
larger initial estimation errors. Similar to iT.jc-filtcring problems, 
we have in general not a bounded positive definite solution of the 
corresponding Riccati differential equation, which does exist for the 
extended Kaiman filler if ilie system satisfies certain observability 
conditions [1]. Tlierefore, a careful choice of the constants in this 
Riccati differential equation is of particular interest. 
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Fractional-Order Systems and FI^D^-ContrQ\\c:rs 
Jgor Podlubny 



.Abstract — Dynamic systems of an arbitrarj' real order (fractional- 
order systems) are considered. A concept of a fractional-order PI^D*'- 
controlJer, involving fractional-order integrator and fractional -order dif- 
ferentiator, is proposed. The Laplace transform formula for a new 
function of the Mittag-Leffler-type made it possible to obtain explicit 
analytical expressions for the unit-step and unit-Impulse response of 
a linear fractional-order system with fractional-order controller botJi 
for the open and closed loop. An example demonstrating the use of 
the obtained formulas and the advantages of the proposed PI^D'^- 
controllers is given. 

Index TfrnM— Fractional-order controllers, fractional-order systems, 
fractional differential equations, Laplace transforms, transfer functions. 



I. Introduction 
Recently, .*;everal authors have con.sidered mechanical systems de- 
scribed by fractional-order state equations [4], [5], [15]. which means 
equatiom involving so-called fractional derivatives and integrals (for 
the introduction to this theor>' see [21]). 
There are also several recent applications in electricity. 
Le Mehaute and Crepy [14] have proposed a concept of a frac- 
to nee — a new electrical circuit element, which has intermediate 
properties between resistance and capacitance. Such a device has been 
experimenlally studied, for example, by Nakagawa and Sorimachi 
[18], A circuit proposed by Oldham and Zoski [22] provides another 
example of a fractance. 

A new capacitor theory developed by Westerlund [28] is ba,scd on 
the use of fractional derivatives. 

New fractional derivative-based models are more adequate than the 
previously used integer-order models. This has been demonstrated, for 
instance, by Caputo [7], Nonnenmacher and Giockle [19], Friedrich 
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